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Abstract Turbulence is ubiquitous in astrophysics. It radi- 
cally changes many astrophysical phenomena, in particular, 
the propagation and acceleration of cosmic rays. We present 
the modern understanding of compressible magnetohydro- 
dynamic (MHD) turbulence, in particular its decomposition 
into Alfven, slow and fast modes, discuss the density struc- 
ture of turbulent subsonic and supersonic media, as well as 
other relevant regimes of astrophysical turbulence. All this 
information is essential for understanding the energetic par- 
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tide acceleration that we discuss further in the review. For 
instance, we show how fast and slow modes accelerate en- 
ergetic particles through the second order Fermi accelera- 
tion, while density fluctuations generate magnetic fields in 
pre-shock regions enabling the first order Fermi accelera- 
tion of high energy cosmic rays. Very importantly, however, 
the first order Fermi cosmic ray acceleration is also possible 
in sites of magnetic reconnection. In the presence of turbu- 
lence this reconnection gets fast and we present numerical 
evidence supporting the predictions of the Lazarian & Vish- 
niac (1999) model of fast reconnection. The efficiency of 
this process suggests that magnetic reconnection can release 
substantial amounts of energy in short periods of time. As 
the particle tracing numerical simulations show that the par- 
ticles can be efficiently accelerated during the reconnection, 
we argue that the process of magnetic reconnection may be 
much more important for particle acceleration than it is cur- 
rently accepted. In particular, we discuss the acceleration 
arising from reconnection as a possible origin of the anoma- 
lous cosmic rays measured by Voyagers as well as the origin 
cosmic ray excess in the direction of Heliotail. 

Keywords Turbulence • Magnetic reconnection • Accelera- 
tion • Cosmic rays 



1 Introduction 

It is well know n that astrophysical fluids are magnetized and 
turbulent (see Armstrong and Wo ol 1 198 It lArmstrong et al 



1995; Verdini and Velli 2007; Lazarian 2009, and references 
therein). For interstellar medium Fig. Q] illustrates the turbu- 
lent power density plotted against the inverse of the scale 
length, with data at large scales, i.e. at small wave numbers 
q expanded using the Wisconsin H a Mapper (WHAM) data 
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amplitude instabilities, from which turbulent motions read- 
ily develop. 

For interstellar medium the drivers of turbulence include 
supernovae explosions that shape the interstella r medium 
(hvIcKee and Ostrikerll977tlNakamura et al.l2 006). accretion 
flows (Kle ssen and Henne belle 2010), magneto-rotational in - 
stability in the galact ic disk ( Sellwood and Balbusl 1999). 



therm al instability (see Kritsuk and Norman 2002; Kov ama and i. Inutsuka 



2002), collimated outflows (see Nakamura and Li 2007), etc.. 
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Fig. 1 Turbulence in the interstellar gas as rev ealed by electron density 
fluctuations. "Big Power Law in the Sky" in lArmstrong et all i ll 9951) 
extended using WHAM data. The slope corresponds to that of Ko l- 
mogorov turbulence. Modified from Chepurno v and Lazarianl (2010). 



Cluster mergers, jets of radio galaxies and turbulent galaxy 
wakes a re among the driver s of turbulence in clusters of 
galaxies (Brun etti et al .120091) . Similarly, the fast plasma flow 
and plasma instabilities provide the natural environment for 
turbulence to develop in the solar wind. In more special cir- 
cumstances, instabilities in the fluid of cosmic rays as well 
as energy injection via magnetic reconnection can serve as 
the source of turbulence. The latter provide special cases that 
we also consider in the paper, but the generic set up is that 
both the evolution of cosmic rays and reconnection is con- 
sidered in the presence of pre-existing turbulence which is 
created by external sources^. 

The model that has been shown to successfully repro- 
duce the pr operties of incompressible MH D turbulence was 
proposed bv lGoldreich and Sridharl d 19951 henceforth GS95) . 



Numerical studies of MHD turbulence in Cho and Lazarian 



on el ectron density fluctuation^] dChepurnov and Lazarianl 
2010h . A similar Kolmogorov- type power law is measure d 

Coleman! d 19681) 



with in situ measurements in Solar wind 



(2002, 2003) allowed to extend the latter model to compress- 
ible fluids. This allowed to change our understanding of the 
effects of turb ulence on the scattering and acceleration of 
cosmic rays (lYan and Lazarianll20 02. 2004). 

Turbulence is known to affect most of the properties of 
fluids, e.g. thermal conductivity, propagation of waves and 



J,™oll2l ^ ii ii r i • i views and books by Moftatt 1978; Dmitruket al. 2001; Schh 

( 1998|j. On much larger scales, clusters of galaxies also re- 
veal turbulent spectrum (Enfi lin and Vo gt 2006). 



_ , Schlickeiser 
Cranmer and van Ballegooijen 2005 



The fact that turbulence is ubiquitous in astrophysical 
environments is the most natural and expected. Magnetized 
astrophysical plasmas generally have very large Reynolds 
numbers due to the large length scales involved, as well as 
the fact that the motions of charged particles in the direction 
perpendicular to magnetic fields are constrained. Laminar 
plasma flows at these high Reynolds numbers R = VLf/v, 
where V and Lt are the velocity and the scale of the flow, 
V is fluid viscosity, are prey to numerous linear and finite- 

1 A more direct evidence comes from the observations of spec- 
tral lines. Ap art from showing non-thermal Doppler broadening (see 
lLarson|[l98lh . they also reveal spectra of supersonic turbulent veloc- 
ity fluctuations when analyzed with techniques like Velocity Channel 
Ana lysis (VCA) of Velocity Coordinate Spect r um (VCS) developed 
(see lLazarian and Pogosvan[ |2000, 2004 1 , 1200(1 12008|) and applied to 
the ob servational data (see lPadoan et alJ i2004, 2009; Chepurno v et alj 

EoU. 

2 More discussion of t h e solar wind turbulen ce can be found in 
iBruno and Carbond j2005T) : lPetrosvan et all d2010l) 



| 2003HVishniac et alT2003l 

lLongairl2010i and references therein). For instance, the cos- 
mic ray streaming instability that is the basic textbook pro- 
cess is being suppres sed in the presence of turbulence in the 

interstellar medium d Yan and Lazarianl2002 : Farmer and Goldreichl 
l2004l[Beresnvak and Laz arian 2008). Similarlv. lLazarian and Vishniad 
(119991 henceforth LV99) shows that the constrains on the 
classical Sweet-Parker reconnection are being lifted in the 
presence of turbulence and the reconnection rate becomes 
fast, i.e. independent on resistivity. 

Fast magn etic reconnection predicted i n the presence of 
turbulence in Lazari an and Vis hniac (1999) e ntails many as- 
trophysical consequences (see e.g. lLazarianluOl U for the 
consequences for star formation). It entails, for instance, the 
efficient accel eration of cosmic rays in reconnection re gions 
as d iscussed in de Go uveia Dal Pino and Lazarianl (12005b (see 
also lLazarianll2005l) . 



3 The relevant analogy: turbulence that is experienced by the pas- 
sengers of an airplane does not arise from the airplane motion, but 
originates due to complex processes taking place in the atmosphere. 
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In what follows, we discuss the nature of MHD turbu- 
lence in §2, its consequences for magnetic reconnection in 
§3 and for cosmic ray acceleration in §4. We describe the 
first order Fermi acceleration that takes place in the recon- 
nection sites in §5. The discussion and summary are pre- 
sented in §6 and §7, respectively. 



2 MHD turbulence in astrophysical plasmas 

MHD turbulence is widely spread process that is essential 
both for astrophysical reconnection and cosmic ray physics. 
The balanced MHD turbulence is relatively well explored 
subject, but other regimes of MHD turbulence, e.g. imbal- 
anced MHD turbulence, viscosity dominated MHD turbu- 
lence are still insufficiently studied and the consequences of 
these regimes are poorly explored. 



2.1 Applicability of MHD approximation 

The paper is focused on the properties of MHD turbulence 
as we claim that the small scale plasma effects are of sec- 
ondary importance for both the dynamics of turbulence at 
sufficiently large scales as well as for reconnection and ac- 
celeration of cosmic rays in the presence of turbulence. 

The applicability of MHD description to plasmas has 
been revisited recently in lEvink et al.1 (1201 ll henceforth ELV 11) 
There three characteristic length-scales were considered: the 
ion gyroradius p,-, the ion mean-free-path length t m fp,i, and 
the scale L of large-scale variation of magnetic and velocity 
fields. Astrophysical plasmas are in many cases "strongly 
collisional" in the sense that £ m f p i <C P;, which is the case 
for the interiors of stars and accretion disks. In such cases, a 
fluid description of the plasma is valid. In the "weakly colli- 
sional" l m fp i 3> Pi,. The ratio 



I, 



mfp.i 



A V A 
InA c ' 



(1) 



follows from th e standard formu la for the Coulomb collision 
frequency (see Fitzpatrickl2008 , Eq. 1 .25). Here A = AnnX^ 
is the plasma parameter, or the number of particles within 
the Debye screening sphere. When astrophysical plasmas 
are very weakly coupled (hot and rarefied), then A is large, 
e.g. of the order of 10 9 or more for the warm component 
of the interstellar medium or solar wind (see Table [TJ- For 
such ratio the expansion over small ion Larmor radius p, 
provides "kinetic MHD equations" which differ from the 
standard MHD by having anisotropic pressure tensor. 

Plasmas that are not strongly collisional further subdi- 
vided into two classes: "collisionless" plasmas for which 
£mfp,i 3> L, the largest scales of interest, and "weakly col- 
lisional" plasmas for which L 3> i m fp,i- I n me latter case 



the "kinetic MHD" description can be further reduced in com- 
plexity at scales greater than £ m fp,i (see ELV11). This re- 
produces a fully hydrodynamic MHD description at those 
scales, with anisotropic transport behavior associated to the 
well-magnetized limit. Among our examples in Table 1 above, 
the warm ionized ISM is "weakly collisional", while post- 
CME current sheets and the solar wind impinging on the 
magnetosphere are close to being "collisionless." 

Additional important simplifications occur if the follow- 
ing assumptions are satisfied: turbulent fluctuations are small 
compared to the mean magnetic field, have length-scales 
parallel to the mean field much larger than perpendicular 
length-scales, and have frequencies low compared to the ion 
cyclotron frequency. These ar e standard assumptions of the 
Goldreich and Sridharl dl995h theory of MHD turbu lence. 



They are the basis of the "gyrokinetic approximation" (Schekochihi n et al 
20071 l2009h . At length-scales larger than the Larmor ra- 
dius p,, another reduction takes place. The incompressible 
shear- Alfven wave modes exhibit dynamics independent of 
compressive motions and can b e described by the "Reduced 
MHD" (RMHD) equations (see|straussll976l:IZank and Matthaeusl 
1992allGoldreich and Sridharl 1995llCho and Lazarianl2003h . 
This fact is essential for the LV99 justifying the use of the 
treatment based on an incompressible MHD model. 



2.2 MHD turbulence: cascades of Alfenic, slow and fast 
modes 

Magnetized turbulence is a tough and complex problem with 
many excellent mo nographs and re views devoted to differ- 
ent aspect of it (see Biskamp 2003, and references therein). 
A broad outlook on the astrophysical im plications of the tur- 
bulenc e can be found in a review by Elmegreen and Scalol 
( 2004), while the effects of turbulence on molecular clouds 
and star formation are well reviewed in McKee and Ostriker 



d2007l) . However, the issues of turbulence spectrum and its 
anisotropies, we feel, are frequently given less attention than 
they deserve. 

We deal with magnetohydrodynamic (MHD) turbulence 
which provides a correct fluid-type description of plasma 
turbulence at large scales. Astrophysical turbulence is a di- 
rect consequence of large scale fluid motions experiencing 
low friction. The Reynolds numbers are typically very large 
in astrophysical flows as the scales are large. As magnetic 
fields decrease the viscosity for the plasma motion perpen- 
dicular to their direction, Re numbers get really astronomi- 
cally large. For instance, Re numbers of 10 10 and larger are 
very common for astrophysical flows. For so large Re the 
inner degrees of fluid motion get excited and a complex pat- 
tern of motion develops. 

While turbulence is an extremely complex chaotic non- 
linear phenomenon, it allows for a r emarkably simple statis- 
tical description (see lBiskampll2003l) . If the injections and 
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Table 1 Parameters for Some Astrophysical Plasmas (adopted from Evin k et alj201lh . 



Parameter 


warm ionized 


post-CME 


solar wind at 




ISM" 


cuiTent sheets^ 


magnetosphere' 


density n, cm ^ 


.5 


7 x 10' 


10 


temperature T , &V 


.7 


10 3 


10 


plasma parameter A 


4x 10 9 


2x 10 10 


5 x 10 10 


ion thermal velocity \>th,u cm/s 


10 6 


3x 10 7 


5x 10 6 


ion mean-free-path £ m fp,i, cm 


6x 10 11 


10'° 


7x 10 12 


magnetic diffusivity A, cm 2 /s 


10 7 


8x 10 2 


6x 10 s 


magnetic field B, G 


io- 6 


1 


io- 4 


plasma beta j8 


14 


3 


1 


Alfven speed v a, cm/s 


3x 10 5 


3x 10 7 


7x 10 6 


ion gyroradius p; , cm 


10 s 


3x 10 3 


6x 10 6 


large-scale velocity U, cm/s 


10 6 


4x 10 6 


5x 10 6 


large length scale L, cm 


10 20 


5x 10'° 


10 8 


Lundquist number Si = ^j- 


3x 10 18 


2x 10 15 


10" 


resistive length* (jj , cm 


5x 10 5 


1 


20 



iNorman and Ferrard dl996t ): lFerrierd i200lh b Bemporad J2008h IZimbardo etZlfcOlOl) 

*This nominal resistive scale is calculated from ijj ~ L(va/U)S l 3//4 , assuming GS95 turbulence holds 

down to that scale, and should not be taken literally when £j: < pj. 



sinks of the energy are correctly identified, we can describe 
turbulence for arbitrary Re and Rm. The simplest descrip- 
tion of the complex spatial variations of any physical vari- 
able, X(r), is related to the amount of change of X between 
points separated by a chosen displacement 1, averaged over 
the entire volume of interest. Usually the result is given in 
terms of the Fourier transform of this average, with the dis- 
placement 1 being replaced by the wavenumber k parallel to 
1 and |k| = 1/|1|. For example, for isotropic turbulence the 
kinetic energy spectrum, E(k)dk, characterizes how much 
energy resides at the interval k, k + dk. At some large scale 
L (i.e., small k), one expects to observe features reflecting 
energy injection. At small scales, energy dissipation should 
be seen. Between these two scales we expect to see a self- 
similar power-law scaling reflecting the process of non-linear 
energy transfer. 

Thus, in spite of its complexity, the turbulent cascade is 
self-similar over its inertial range. The physical variables are 
proportional to simple powers of the eddy sizes over a large 
range of sizes, leading to scaling laws expressing the depen- 
dence of certain non-dimensional combinations of physical 
variables on the eddy size. Robust scaling relations can pre- 
dict turbulent properties on the whole range of scales, in- 
cluding those that no large-scale numerical simulation can 
hope to resolve. These scaling relations are extremely im- 
portant for obtaining an insight of the processes at the small 
scales. 

The pre sence of a magnetic field makes MHD turbulence 
anisotropic (Montgomery and Turner! 1 98 ltuMatthaeus et al.1 



1983 



1995 



Sheb alin et al.1 1983l:lHigdonTl 984 



Goldreich and Sridhar 



see lOughton et al.L 120031 for a review). The relative 
importance of hydrodynamic and magnetic forces changes 
with scale, so the anisotropy of MHD turbulence does too. 



Many astrophysical results, e.g. the dynamics of dust, scat- 
tering and acceleration of energetic particles, thermal con- 
duction, can be obtained if the turbulence spectrum and its 
anisotropy are known. As we discuss below, additional im- 
portant insight can be obtained if we know turbulence inter- 
mittency. 

Estimates of turbulence anisotropy obtained in relation 
to the observations of magneti c fluctuation of the outer he- 
liosphere and solar wind (see IZank an d Mattha euslll992b . 
and references therein) provided, for an extended period of 
time, the only guidance for theoretical advances. This re- 
sulted in a picture of MHD turbulence consisting of 2D "re- 
duced MHD" perturbations carrying approximately 80% of 
energy and "slab" Alfvenic waves carrying the remaining 
20% of energy (see iMatthaeus et al. 20021 and references 
therein). In other words, in the suggested picture the MHD 
turbulence was presented by two anisotropic components, 
one having wave vectors mostly perpendicular to magnetic 
field (the 2D one), the other having them mostly parallel 
to magnetic field (the slab one). This model became a de- 
fault one for many calculations of the propagation of cosmic 
rays (see lBieber et al.lll9 88. 1994). On the contrary, guided 
mostly by compressible MHD numerical simulations, the 
interstellar community adopted a model of the MHD tur- 
bulence where the basic MHD modes, i.e. slow, fast and 
Alfvenic are well coupled together and efficiently di ssipate 
energy in shocks dstone et al.ll 19981 Mac Lowlll999l) . Little 
cross-talk between the two communities did not stimulate 
the interdisciplinary debates on the nature of MHD turbu- 
lence, which was regretful, as the heliospheric community 
has the advantage of the in-situ spacecraft measurements. 

In spite of the intrinsic limitations of the "brute force" 
approach, we feel that reliable results can be obtained nu- 
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merically if the studies are focused on a particular property 
of turbulence in order to get a clear picture of the underly- 
ing physics occurring on small scales ("microphysics") that 
cannot be resolved in "global" interstellar simulation^]- 



Fo r instance, numerical studies in lCho and Lazarianl (2002. 
2003b showed that the Alfvenic turbulence develops an in- 
dependent cascade which is marginally affected by the fluid 
compressibility. This observation corresponds to theoretical 
expectations of the GS95 theory that we brie fly describe be- 
low (see also Lithwick and GoldreichlEoOl ). In this respect 
we note that the MHD approximation is widely used to de- 
scribe the actual magnetized plasma turbulence over scales 
that are much larger than both the mean free path of the par- 
ticles and their Larmor radius (see Kulsrudfl983 , 2005 , and 
references therein). More generally, the most important in- 
compressible Alfvenic part of the plasma motions can de- 
scribed by MHD even below the mean free path but on the 
scales larger than the Larmor radius. 

We claim that, while having a long history of ideas, the 
theory of MHD turbulence has become testa ble recently due 
to the advent of num erical simulations (see Biskamp 2003) 
which confirmed (see Cho and Lazarianl2005l and references 
therein) the prediction of magnetized Alfvenic eddies being 
elongated in the dire ction of magnetic field (see lShebalin et al.1 
1983: lHigdonll984h and provided results consistent with the 
quantitative relations for the degree of eddy elongation ob- 
tained in GS95. Indeed, GS95 made predictions regarding 
relative motions parallel and perpendicular to B for Alfvenic 
turbulence. The model did not predict the generation of any 
"slab" modes and, instead of pure 2D Alfvenic modes, pre- 
dicted that most of the Alfvenic energy is concentrated in 
the modes with a so-called "critical balance" between the 
parallel and perpendicular motions. The latter can be under- 
stood within intuitive picture where eddies mixing magnetic 
field perpendicular to its loca|^| direction induce Alfvenic 
waves with the period equal to the period of the eddy ro- 
tation. This results in the scale-dependent anisotropy of ve- 
locity and magnetic perturbations, with the anisotropy being 
larger for smaller eddies. 

The hydrodynam ic counterpart of the MHD turbulence 
theory is the famous iKolmogorovl (119411) theory of turbu- 
lence. In the latter theory energy is injected at large scales, 
creating large eddies which correspond to large Re num- 
bers and therefore do not dissipate energy through viscos- 

4 By contrast, numerical simulations that deal with many physical 
conditions simultaneously cannot distinguish between the effects of 
different processes. Moreover, they inevitably have a more restricted 
interval of scales on which energy is injected by numerics, initial con- 
ditions, or boundary conditions. Their results are, therefore, difficult to 
interpret in physical terms. 

5 The notion of the direction being local is critical. Small eddies are 
affected by magnetic field in their vicinity, rather than a global field. 
No universal scalings are possible to obtain in the frame of the mean 
magnetic field. 



ity^but transfer energy to smaller eddies. The process con- 
tinues until the cascade reaches the scales which are small 
enough that the energy is dissipated during one turnover 
time of the eddies at the corresponding scales. In the ab- 
sence of compressibility the hydrodynamic cascade of en- 
ergy is ~ vf/ tease l — const, where V/ is the velocity at the 
scale I and the cascading time for the eddies of size I is 
t C ask,l ~ l/ v l- From this the well known relation v/ ~ Z 1 / 3 
follows. 

A frequent mental picture that astrophysicists have of 
the Alfvenic turbulence is based of Alfven waves with wave 
vectors along the magnetic field. This is not true for the 
strong Alfvenic turbulence which, similar to its hydrody- 
namic counterpart, can be described in terms of eddieo 
However, contrary to Kolmogorov turbulence, in the pres- 
ence of dynamically important magnetic field eddies become 
anisotropic. At the same time, one can imagine eddies mix- 
ing magnetic field lines perpendicular to the direction of 
magnetic field. 

The nature of Alfvenic cascade is expressed through the 
critical balance condition in GS95 model of strong turbu- 
lence, namely, 

/,7 h 



(2) 



where v/ is the eddy velocity, while the and l± are, re- 
spectively, eddy scales parallel and perpendicular to the lo- 
cal direction of magnetic field. The critical balance condi- 
tion states that the parallel size of an eddy is determined by 
the distance Alfvenic perturbation can propagate during the 
eddy turnover. The notion of local is important, as no uni- 
versal relations exist if eddies are t reated with respect to the 
global mean magnetic field ( LV99; Cho and Vishniac 2000; 



Maron and Goldrei chll200ll: iLithwick and GoldreichlbOQl 
Cho and Lazarianll2002 ). 



The critical balance is the feature of the strong turbu- 
lence, which is the case when the turbulent energy is in- 
jected at Va- If the energy is injected at velocities lower 
than Va the cascade is weak w ith / 1 of the eddies increas- 
ing while l \\ staying the same dNg and Bhattacharieelll996t 
Lazarian and Vishniadl999 ; Galtier et al.l2002h . In other words, 
as a result of the weak cascade the eddies get thinner, but 

6 Reynolds number Re = LfV/v = (V /L f )/{v /L 2 f ) which is the ra- 
tio of an eddy turnover rate Z^ di = V/Lf and the viscous dissipation 
rate zj} — T\/Iy. Therefore large Re correspond to negligible viscous 
dissipation of large eddies over the cascading time z casc which is equal 
to Z e ddy in Kolmogorov turbulence. 

7 The description in terms of interacting wave packets or modes is 
also possible with the corresponding wave vectors tending to get more 
and more perpendicular to the magnetic field as the cascade develops. 

8 To stress the difference between local and global systems here we 
do not use the language of ^-vectors. Wave vectors parallel and per- 
pendicular to magnetic fields can be used, if only the wave vectors are 
understood in terms of a wavelet transform defined with the local ref- 
erence system rather than ordinary Fourier transform defined with the 
mean field system. 
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Table 2 Regimes and ranges of MHD turbulence. 



Type 


Injection 


Range 


Motion 


Ways 


of MHD turbulence 


velocity 


of scales 


type 


of study 


Weak 


v L <v A 


[L, I trans] 


wave-like 


analytical 


Strong 










subAlfvenic 


V L <V A 


[hrans: huin\ 


eddy-like 


numerical 


Strong 










superAlfvenic 


V l >V a 


) lmin\ 


eddy-like 


numerical 



L and /„„„ are injection and dissipation scales 

hrum and Ia are given by Eq. j5) and Eq. J6), respectively. 



preserve the same length along the local magnetic field. This 
decreases l± and eventually makes Eq. (0 satisfied. If the 
injection velocity is Vl and turbulent injection scale is L, 
the transition to the strong MHD turbulence happens at the 
scale 1(vi/Va) 2 an d the velocity a t this scale is V str ong = 
V A (vi/V A ) 2 (LV99: lLazarianll2006l) . Thus the weak turbu- 
lence has a limited, i.e. [/,/(v/ /Va) 2 ] inertial interval and get 
strong at smaller scales. 

While GS95 assumed that the turbulent energy is in- 
jected at Va at the injection scale /, LV99 provided gen- 
eral relations for the turbulent scaling at small scales for the 
case that the injection velocity Vl is less than or equal to V A , 
which can be written in terms of h and l±: 



In Rii 



2/3 



1/3 



M 4/3 

Vl) 



(3) 



1/3 



It is important to stress that the scales l± and /|| are mea- 
sured with respect to the system of reference related to the 
direction of the local magnetic field "seen" by the eddy. 
This notion was not present in the original formulation of 
the GS95 theory and was added to it in LV99. The local 
system of reference was later used in numerical studi es in 



Cho and Vishniac (2000), Maron and Goldreich (2001), and 



Cho et alJ (hool " testing GS95 theory. In terms of mixing 



motions, it is rather obvious that the free Kolmogorov-type 
mixing is possible only with respect to the local magnetic 
field of the eddy rather than the mean magnetic field of the 
flow. 

The original GS95 picture deals with transAlfvenic tur- 
bulence, i.e. with Va ~ Vl- For low Alfvenic Mach numbers, 
i.e. for Va 3> Vi at large scales ~ L the turbulence is weak 



(see Ng and Bhattachariee 1997; 


Lazarian and Vishniac 1999 


Galtier et al.l2000H and magnetic fields are slightly perturbed 



by propagating Alfven waves. The wave packets in weak tur- 
bulence evolve changing their perpendicular scale l±, while 
their scale /|| along the magnetic field does not change. As 
scaling of weak turbulence predicts V\ ~ Vl(1± /L) 1 / 2 (LV99), 
at the scale 



In 



L(V L /V A ) 2 = LM 2 



the critical balance condition 1\\/Va ~ l±/Vi is getting sat- 
isfied making turbulence strong. It is easy to see that the 
velocity corresponding to l tmm is V trans ~ V L {V L /V A ). 

For superAlfvenic turbulence the situation is somewhat 
different. Magnetic field gets dynamically important as soon 
as its energy density exceeds the energy of eddies at the 
Ohmic dissipation scale, which translates into Alfvenic ve- 
locity getting larger than the velocity of eddies at the Ohmic 
dissipation scale or the ion Larmor radius, whichever is larger. 
As the velocity in Kolmogorov turbulence scale as v; ~ Z 1 / 3 , 
it is clear that even weak magnetic field can make a sig- 
nificant impact on the dynamics of the smallest eddies. In 
view of that it is advantageous to introduce a scale at which 
the magnetic field gets dynamically important and the nature 
of the turbulen ce changes from hydrodynamic to MHD (see 
LazarianlboOrjl) . namely, 



(4) l A =L{V A /V L ) 3 =LM j 



-3 



(6) 



The relations predict ed in GS95 were confirmed numeri- 
cally for incompressib le ( Cho and Vishniad200C ; 



2001; Cho etal 



Beresnvak2011 



2003 



20021; iBeresnvak and Lazariarl2009 



) and compressi ble MHD turbulenc 



Maron and Goldreich 
20l d; 



2002, 2003 : Kowal and Lazarian 20ld see also Cho et al 



(Cho and Lazarian 



for a review). They are in good agreement with ob- 
served and inferred astrop hysical spectra. A remarkable fact 
revealed in lCho et al.1 (120021) is that fluid motions perpendic- 
ular to B are identical to hydrodynamic motions. This pro- 
vides an essential physical insight and explains why in some 
respects MHD turbulence and hydrodynamic turbulence are 
similar, while in other respects they are different. 

GS95 provided theoretical arguments in favor of weak 
coupling between fast and Alfven mo des, and low impact of 
slow m odes to Alfven modes (see also lLithwick and Goldreich 
200 ll) . This challenged the paradigm accepted by the inter- 
stellar community. While the decomposition of MHD per- 
turbations into fundamental MHD waves was widely used in 



(5) 



9 So me studies of MHD compressible turbulence, e.g. Vestu to et alj 
(2003) did not perform a decomposition of MHD perturbations into 
Alfve n, slow and fast modes as it is done in lCho and Lazariar] i2002L 
2003). They did not use local system of reference for which the GS95 
scaling is fonnulated. Therefore a direct comparison of their results 
with the GS95 predictions is difficult. 
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the literature (see Dobrowolnv et al . 1980) it was usually as- 
sumed that the Alfvenic waves exist and int eract with othe r 
waves for many periods (see a discussion in IZweibellEooll) . 
This is not the case of the GS95 model of turbulence, where 
the Alfven modes non-linearly decay within one wave pe- 
riod. This reduces the time of interaction and therefore the 
coupling. Interestingly enough, in GS95 model, the Alfvenic 
modes can affect slow modes, but the opp osite is not true. 
These result s we re successfully tested in 



Cho and Lazarian 



§2.1), one should keep in mind that the strength of Alfvenic 
interactions increases with the decrease of the scale along 
the cascade. Therefore the Alfvenic turbulence gets eventu- 
ally strong, while both the inertial range and the astrophys- 
ical utility of the weak Alfvenic cascade are limited. The 
interaction of weak Alfvenic turbulence with fast modes has 
dependences on the angle between B and the wave vector 
dChandranll2005l) . 

As the Mach number of turbulence increases, shocks 
should play more important role in the dynamics of com- 
pressible motions. The transition from the regime when com- 
pressible motions can be well described by the fast and slow 
ode cascades versus the situation when shocks carry an ap- 



Beresnvak et alj 02OO5L henceforth BLC05) shows that in 



J2002L l2003h and lKowal and Lazarianl (bold 

Some of the relevant results are illustrated in Figure [2] 
Contrary to the common expectation, the modes exhibited 
nice scaling laws that allow further analytical and numerical rn 

applications For instance, numerical studies inlch o and Lazaria j^r eciable part of the energ y deserves more studies. A study 
J20021 l2003h revealed that the GS95 scalings are valid for 
the Alfvenic part of the turbul ence cascade even in the highly 
compressible regime (see also Beresnvak and Lazarianl2006 ). 

The statistical decomposition of MH D turbulence into 
Alfven , slow and fast modes suggested in lCho and Lazarian 
(2002) was successfully test ed for slow modes of ma gnet- 
ically dominated plasmas in Cho and Lazarianl ( 2003 ). The 
decomposition of MHD turbulent motions into the discussed 
thr ee cascades of fundamental modes have been confirmed 
by iKowal and Laza rian l l2010l) . where in order to improve 
the decomposition, the wavelet transformation were used. 
The large scale components (large with respect to the size 
of the local wavelet), which determine the vector base of 
mode projection in Fourier space, are obtained by averaging 
over the size of the wavelet. Therefore, the wavelet trans- 
formations allow to approach closer to the decomposition 
in the local system of reference than the earlier statistical 
approach based on Fourier analysis could do. In Figure [3] 
we show that the new wavelet decomposition results (upper 
row) are very similar to the ones obtained with the Fourier 
method (lower row). At the same time, the study of spectra 
and anisotropy of velocity, magnetic field reveals the advan- 
tages of the wavelets in comparison with the Fourier tech- 
nique for studying turbulence with weak mean field. For tur- 
bulence with B mean ~ SB the obtained r esults are consistent 
with the Cho and Lazarianl ( 2002 , 2003 ) studies. 

Nevertheless, it would be wrong to say that we have a 
complete understanding of the scaling of MHD modes and 
their interactions. First of all, one should distinguish weak 
and strong Alfvenic turbulence. The weak turbulence is es- 
sentially 2E0 , with the turbulent cascade creating more 
structure perp endicular to magne tic field as the turbulence 
cascades (see Galtier et al. 2000h . Such a cascade emerges 
when the driving of turbulence at the outer scale is weak, 
i.e. the injection velocity is much less than the Alfven veloc- 
ity. Although the weak turbulence picture corresponds to the 
early representation of MHD turbulence (see discussion in 



Practica l studies of non-linearity of tu r bulence in Solar wi nd are 
discussed in Sahr aoui and Go ldstein (2010); Ghosh et al J 1 201 1 ). 



the space between shocks the description of the perturba- 
tions with slow modes is valid for all the Mach numbers 
studied (up to 10). 

The total velocity field contains two components: solenoidal, 
which is equivalent to the incompressible part, and potential, 
which contains the compressible part of the field and the re- 
maining part which is curl and divergence free. In Table [3] 
we show the percentage contribution of each component to 
the total velocity field. We see that the compressible part 
constitutes only a fraction of the total field. However, the 
magnitude of this fraction is different for sub and supersonic 
models. In the case of sub- Alfvenic turbulence, it is about 
3% in subsonic models and about 7% in supersonic mod- 
els, which confirms a higher efficiency of the compression 
in the presence of supersonic flows. Furthermore, the frac- 
tion also changes when we compare models with strong and 
weak magnetic fields. The velocity field, in the presence of 
a weak magnetic field, contains about 5% of the compress- 
ible part in the model with j$ s ~ 0.7 and even up to 16% in 
models with j$ s > 1 . The consequence of the presence of a 
strong magnetic field results in a reduction of the compress- 
ible part of the velocity field by a factor of 2. This indicates 
a substantial role of the magnetic field in the damping of the 
generation of the compressible flows. 

The wavelet decomposition, important for MHD turbu- 
lence, separates the velocity field into three different MHD 
waves: an incompressible Alfven wave and slow and fast 
magneto acoustic waves, of both which are compressible. In 
Table [3] we included the percentage amount of these com- 
ponents in the total velocity field. As we see, most of the 
energy is contained in the Alfven wave. It is almost 60% in 
the case of sub-Alfven turbulence, and about 50% for super- 
Alfvenic turbulence. The slow wave contains approximately 
1/3 of the total energy. However, for the super- Alfvenic case, 
this amount is slightly higher. Table [3] suggests that the slow 
wave is weaker when the turbulence becomes supersonic. 
We do not see a similar behavior for the Alfven wave in 
the case of models with a strong magnetic field. This effect 
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Fig. 2 Statistics of MHD Turbulence: (a). Scaling of compressible motions for plasma with magnetic pressure ten times the gas pressure. This 
regime is important for molecular clouds (from CL02). The velocity (solid line) and B (dashed line) spectra are plotted against k (= l/(eddy size)). 
They show well defined statistical properties that allow further fruitful applications, (b). The simulations of decaying turbulence show that the rate 
of the decay of the total energy is a strong function of the imbalance of the energy contained in waves ini tially moving in opposite directions. 
In the lowest curve, the waves have the same amplitude. The energy decays more slowly with imbalance ICho et aL 2002). (c). The evolution 
of kinetic energy in 3-D compressible turbulence, initially started with V • v = as if the gas were incompressible. The dashed curve shows the 
evolution of the V • v = motions. We see that Alfvenic turbulence creates only a marginal amount of compressible motions, suggesting that 
Alfvenic modes should evolve independently of the compressible cascade (CL02). (d). Magn etic fluctuation s persist beyond the turbulent damping 
scale at large k, while hydrodynamic fluctuations damp out in partially ionized gas {Lazarian et al. 2004). This viscosity-dominated regime of 
turbulence may dominate small scale structure of partially ionized gas. (e-f). Isocontours of equal correlation for Alfven and fast modes (CL02). 
(e) The Alfvenic motions are much more correlated along B than perpendicular to it. (f) In contrast, fast magnetosonic fluctuations show essentially 
circular (isotropic) isocontours of correlation. 



could also take place in the super- Alfvenic models, but it is 
weakened by relatively large errors. An interesting depen- 
dence is observed in the case of the fast wave. Although the 
fast wave is the weakest among all MHD waves, it strongly 
depends on the regime of turbulence. Similarly to the com- 
pressible part of the velocity field, it is stronger for models 
with a weak magnetic field. In addition, it is much stronger 
when turbulence is supersonic, but this strength seems to be 
weakly dependent on the sonic Mach number. 



2.3 Density structure of MHD turbulent flows 

Density structure of turbulence is usually discussed in rela- 
tion to star formation studies. However, we shall show in S|6] 
that fluctuations of density play important role for energetic 



particle acceleration. Thus it is important to survey the basic 
properties of the turbulent density field. 

The power spectrum of density fluctuations is an im- 
portant property of a compressible flow. In some cases, the 
spectrum of density can be derived analytically. For nearly 
incompressible turbulent motions in the presence of a strong 
magnetic field, the spectrum of density scales similarly to 
the pressure, i.e. E p (k) ~ kr 1 ^ if we conside r the poly- 
tropic equation of state p = apY ( Biskamp1l2003 ). In weakly 
magnetized nearly incompressible MHD turbulence, how- 
ever, velocities convect density flu ctuations passively induc - 
ing the spectrum E p (k) ~ k~ 5 ^ (Montgomery et al.l ll987). 
In supersonic flows, these relations are not valid anymore 
because of shocks accumulating matter into the local and 
highly dense structures. Due to the high contrast of density, 
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Fig. 3 The topology of the Alfven, slow and fast mode (left, middle, and rig ht column, respectively) i n the plane p arallel to the mean field . 
The upper and lower rows show the MHD modes separated using the wavelet iKowal and LazarianlFioioh and Fourier iCho and Lazariarj|2003h 
methods, respectively. Texture shows the direction of velocity field multiplied by the local amplitude. 



Table 3 Percentage Amount of the Kinetic Energy Contained Within Each Velocity Component in Kowal and Lazarian (2010). Errors correspond 
to a measure of the time variation. 





^lncomp. 


^comp. 


V A 


Vs 


V f 


-0.8 


-0.7 


96.5 ±() « 


3.3 ±tvli 


58 ±4 


37 ±J 


48 ±o.y 


-2.5 


-0.6 


93 ±2 


7 ±2 


58± 5 


33 ±4 


9 ±2 


-7.5 


-0.5 


92 ±2 


7 ±2 


56± 4 


36 ±4 


8.0 ±( " 


-0.7 


-2.1 


95 ±2 


5± 2 


52± 4 


42 ±4 


6.2 ±as 


-2.5 


-2.1 


86±> 


1 4 ±2 


47 ±3 


37 ±4 


16 ±2 


-7.5 


- 1.9 


84 ±2 


16 ±2 


47 ±4 


33 ±4 


20 ±2 



the linear relation Sp = c 2 s 8p is no longer valid, and the 
spectrum of density cannot be related to pressure straight- 
forwardly. 

In Figure|4]we present the power spectra of fluctuations 
of density for models with different ^# s . As expected, we 
note a strong growth of the amplitude of density fluctuations 
with the sonic Mach number at all scales. This behavior is 
observed both in sub-Alfvenic as well as in super- Alfvenic 
turbulence (see Fig.|4|. In Table|4]the spectral index of den- 
sity a p and the logarithm of density a\ ogp is presented. The 
width of the inertial range is shown by the range of solid 
lines with slopes —5/3 and —2 in all spectra plots. Table [4] 
also shows the errors of estimation which combine the er- 
ror of the fitting of the spectral index at each time snapshot 



and the standard deviation of variance of 0Cp.i O gp in time. 
The slopes of the density spectra do not change significantly 
with ./# s for subsonic experiments and correspond to ana- 
lytical estimations (about —2.2, which is slightly less than 
—7/3, for turbulence with j$a ~ 0.7 and about — 1 .7, which 
is slightly more than —5/3, for weakly magnetized turbu- 
lence with ~ 7.3). Such an agreement confirms the va- 
lidity of the theoretical approximations. Those, nevertheless, 
do not cover the entire parameter space. While the fluid mo- 
tions become supersonic, they strongly influence the den- 
sity structure, making the small-scale structures more pro- 
nounced, which implies flattening of the spectra of density 
fluctuations (see value s for j$ s > 1.0 in Table [U see also 
Beresnvak et aDl2005h . 
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Fig. 4 Spectra of density for experiments with different values of and with ^4 ~ 0.7 (left panel) and ~ 7.3 (right panel) for models with 
medium resolution (256 3 ). Grey error bars signify the varian ce of spectra with t ime. The solid lines with slopes —5/3 and —2 cover the inertial 
range used to estimate the spectral indices of p spectra. From lKowal et alj (2007). 

Table 4 Slopes of the Power Spectrum of Density and the Logarithm of Density Fluctuations in Kowa feTal] 120071) . 





je A ~ 0.7 






.^A~7 




J( s 


a p 


^logp 


Jt s 


dp 


^logp 


023 ±o.oi 


-2.3 ±a3 


-2.3 ±0 - 3 


0.26 ±003 


_ 17 ±0.3 


_ 17 ±0.3 


0.33± 001 


-2.2 ±a3 


_2 2 ±0 - 3 


0.36± 004 


_ 17 ±0.3 


_ 17 ±0.3 


0.68 ±003 


-2.0 ±0 - 3 


_ 21 ±o.3 


074 ±0.06 


-1.6 ±0 - 2 


-1.6 ±0 - 3 


220 ±o.03 


-1.3 ±a2 


-2.0 ±() ' 2 


2, 34 ±0.08 


_ 12 ±0.2 


-1.6 ±a2 


7.0 ±a3 


-0.5 ±al 


_ 17 ±o.2 


71 ±0.3 


-0.6 ±a2 


-1.5 ±a2 



The values have been estimated within the inertial range for models with ~ 0.7 (left) and jMa ~ 7 (fight). Errors of spectral indices combine 
the errors of estimation at each time snapshot and the standard deviation of variance in time. Errors for sonic Mach numbers are the standard 
deviation of their variance in time calculated over the period starting from t > 5 to the last available snapshot. 



Another question is the anisotropy of density and the 
logarithm of density structures. For subsonic turbulence it is 
natural to assume that the density anisotropics will mimic 
velocity a nisotropie s in th e GS95 picture. This was con- 
firmed in Icho et al] d2003lh who observed that for super- 
sonic turbulence the contours of densit y isocorrelation be- 



comes round, corresponding to isotropv. lBeresnvak et al.l (|2005) 
however, showed that anisotropies restore the GS95 form if 
instead of density one studies the logarithm of density. This 
is due to the suppression of the influence of the high density 
peaks, which arise from shocks. It is these peaks that mask 
the anisotropy of weaker, but more widely spread density 
fluctuations. 

In Figure [5] we show lines that mark the correspond- 
ing separation lengths for the second-order structure func- 
tions parallel and perpendicular to the local mean magnetic 
fieldFhn the case of subAlfvenic turbulence, the degree of 

1 1 The local mean magnetic field was computed using the procedure 
of smoothing by a 3D Gaussian profile with the width equal to the sep- 
aration length. Because the volume of smoothing grows with the sep- 
aration length /, the direction of the local mean magnetic field might 
change with / at a n arbitrary point. This is an extension of the proce- 
dures employed in lCho et alj 120020 . 



anisotropy for density is very difficult to estimate due to the 
high dispersion of points. However, rough estimates suggest 
more isotropic density structures, because the points extend 
along the line l« ~ For models with j^a ~ 7.3 the points 

in Figure[5]have lower dispersion, and the anisotropy is more 

2/3 

like the type from GS95, i.e. /|| ~ / ± . In both Alfvenic 
, regimes, the anisotropy of density does not change signif- 
icantly with ^# s . Plots for the logarithm of density show 
more smooth relations between parallel and perpendicular 
structure functions. The dispersion of points is very small. 
Moreover, we note the change of anisotropy with the scale. 
Lower values of structure functions correspond to lower val- 
ues of the separation length (small-scale structures), so we 
might note that the logarithm of density structures are more 
isotropic than the GS95 model at small scales, but the anisotropy 
grows a bit larger than the GS95 prediction at larger scales. 
This difference is somewhat larger in the case of models 
with stronger external magnetic field (compare plots in the 
left and right columns of Figure[5), which may signify their 
dependence on the strength of B e xt- The anisotropy of logp 
structures is marginally dependent on the sonic Mach num- 
ber, similar to the density structures. All these observations 
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Fig. 5 Left: Anisotropy for the 2'"' -order SF for the logarit hm of density for ^4 ~ 7.3; Right: the same is for .^4 ~ 7.3. The logarithm of density 
traces the anisotropies of turbulent magnetic field. From Kowal et al. (2003)- 



allow us to confirm the previous studies (see lBeresnvak et al. 

that suggest that the anisotropy depends not only on 
the scale but on ^a- 

Probability distribution functions (PDFs) of density and 
column density give us information about the fraction of 
the total volume occupied by a given value of a measured 
quantity. For the case of compressible turbulence in which 
evolution is described by the Navier-Stokes equation and 
the isothermal equa tion of state, the PDF of density ob eys 
a lognormal form ( Passot and Vazquez-Semadenilll998 ). If 
we select a point in space and assume that the density at this 
point results from subsequent events perturbing the previ- 
ous density, then the final density is a product poi7,(l + 5/), 
where po is the initial density at the selected point and 5,- 
is a small compression/rarefaction factor. By the power of 
the Central Limit Theorem, the logarithm of the resulting 
density, log(p / po) = log(l + 5,-), should obey a Gaussian 
distribution. 



It was discussed in Vazquez-Semadeni and Garcfal ( 2001 ) 
that a parameter, which we for our convenience denote N , 
can determine the form of the column density PDF in molec- 
ular clouds. This parameter is the ratio of the cloud size to 
the decorrelation length of the density field. They defined 
the decorrelation length as the lag at which the density auto- 
correlation function (ACF) has decayed to its 10% level. If 
the density perturbation events are uncorrelated for 8 > 1, 
large values of the ratio K imply that the Central Limit The- 
orem can be applied to those events. 

The decorrelation length estimated from the ACF of the 
density for our models ranges from about 20 cell sizes of the 
computational mesh for supersonic models to about 50 cell 
sizes for subsonic models, which corresponds to 5 < K < 13 
if we take the size of the computational box as the column 
length (in the case of medium resolution, it is equal to 256 
cell sizes). Found values of K signify that in our models at 
least partial convergence to a Gaussian PDF should occur. 



In the top row of Figure [6] we show PDFs of density 
normalized by its mean value for all models with medium 
resolution (256 3 ). The plot on the left shows results obtained 
from sub-Alfvenic experiments (^Ma ~ 0.7), and the plot on 
the right results from super-Alfvenic models (./^a ~ 7.3). 

Our plots confirm the strong dependence of PDFs of 
density on the sonic Mach number, an already known and 

well-studied proper ty of density fluctuations in compress- 

ible tu rbulence (see lVazquez-Semadeni and GarcfabOQll Ostriker 



2003, and references therein). For most of the models PDFs 



are lognormal functions. For super-Alfvenic turbulence PDFs 
are very symmetric about a vertical line crossing their max- 
ima. However, for models with a strong external magnetic 
field (sub-Alfvenic turbulence, ^a ~ 0.7) and very low pres- 
sure (supersonic turbulence, > 1.0), the shape of the 
density PDFs is significantly deformed, and its lower value 
arm ends in higher densities than in the case with a weak 
magnetic field (compare models for Jt s ~ 7 drawn with 
solid lines in the left and right plots of Figure |6j. This sup- 
ports a hypothesis that the gauge symmetry for logp that 
exists in Navier-Stokes equations is broke n in MHD equa- 
tions because of the magnetic tension term ( Beresnva k et alJ 



2005), which physically manifests itself by preventing the 



formation of highly underdense regions. In the higher den- 
sity part of the distribution we do not see a similar effect, 
because the highly dense structures are created mainly due 
to shocks, and in this case, the magnetic tension is too weak 
to prevent the condensation. 

In Figure [6] we included the degree of variation of PDFs 
in time as grey error bars. We note that the departure of 
PDFs from their mean profiles is very small in the middle 
part around the mean value. The strongest time changes are 
observed in the low- and high-density tails, but the PDFs for 
different models are still separable. 

According to lCho and Lazarian d2003h . the relations be- 
tween the variance of density fluctuations and the sonic Mach 
number are 8p /po ~ when the magnetic pressure dom- 
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Fig. 6 Normalized probability distribution functions (PDFs) of density p/p (top row) for models with different values of ,/# s and for ~ 0.7 
(left column) and .4t\ ~ 7.3 (right column). Grey error bars signify the standard deviation of the PDFs, show ing the strength of the departure at a 
single time snapshot from its mean profile averaged over time from t > 5 to the last available snapshot. From Kowal et al. (20071). 



inates and 5p/po ~ ^# s 2 when the gas pressure dominates. 
The former case is observed when j$ s ^> j^a- Indeed, in 
FigureQwe see that the mean standard deviation of the den- 
sity fluctuations (8p /po) scales with ./# s almost linearly 
when .4K S > .Ma- When ./#,, < j$a> the relation is much 
steeper. A similar behavior is observed in the case of super- 
Alfvenic turbulence (Fig. [7] right) although the gas pres- 
sure dominates in all models. For models with a very small 
value of the sonic Mach number, the relation is different. 
The value of (8p /po) depends much less on j% s . 



3 Outstanding issues of MHD turbulence theory 

3.1 The actual spectral slope of MHD turbulence 

Recently there has been a debate concerning an asymptotic 
scaling of MHD turbulence. While the GS95 theory pre- 
dicts Kolmogorov's -5/3 sca ling, the modification of the- 
ory in Boldvrev(2005, 2006) claim^l the power law index 
of -3/2. iBeresnvak and Lazarianl (120091 120101) showed that 
MHD turbulence is less local then its hydrodynamic coun- 
terpart and therefore the shallow spectrum in simulations is 
expected as a result of an extended bottleneck intrinsic to 
MHD turbulence. Below we discuss more recent numerical 
tests related to this controversy. 

One of the way to discriminate different models is by 
means of direct numerical simulations. A quantitative method 
called a resolution study is a traditional way to establish a 
correspondance between numerics and theory. This method 
has been used in the studies of hydrodynamic turbulence, 



The model adds to GS95 the effect of the alignment of velocity 
and magnetic field, e.g. lTineetaljft986h and which has been first nu- 
merically quantified in 3D numerical studies IBeresnvak and Lazarianl 
(2006). It requires that the dynamical alignment should increase with 
the decrease of the scale, whic h, however, contradicts to numerical cal- 
culations in Beresnyak (2012). 



e.g. lYeung and Zhoul dl997l) : iGotoh et all (120021) . etc. We 
perform several simulations with different Reynolds num- 
bers. If we believe that turbulence is universal and the scale 
separation between forcing scale and dissipation scale is large 
enough, the properties of small scales should not depend on 
how turbulence was driven. This is because MHD equations 
and their numerical approximations contain no scale explic- 
itly, save for the grid scale, so simulation with a smaller dis- 
sipation scale could be considered, as a simulation with the 
same dissipation scale, which is a multiple of grid scale, but 
larger driving scale. E.g., the small scale statistics in a 1024 3 
simulation should look the same as small-scale statistics in 
512 3 , if the physical size of the elementary cell is the same 
and the dissipation scale is the same. 

The above scaling argument requires that the geometry 
of the elementary cells is the same and the actual numerical 
scheme used to solve the equations is the same. The scal- 
ing argument does not require high precision on the dissipa- 
tion scale or a particular form of dissipation, either explicit 
or numerical. This is because we only need the small-scale 
statistics to be similar in the particular series of simulations. 
This is achieved by applying the same numerical procedure. 
On the other hand, if turbulence is local, this ensures that 
small scales are only marginally affected by large scales. In 
practice, the scaling argument, which is also known as a res- 
olution study, is done in the following way: the averaged 
spectra in two simulations are expressed in dimensionless 
units corresponding to the expected scaling, for example a 
E(k)k 5 ^e~ 2 ^ is used for hydrodynamics, and plotted ver- 
sus dimensionless wavenumber kr\, where dissipation scale 
77 correspond to the same model, e.g. 77 = (v 3 /e) 1 / 4 is used 
for scalar second order viscosity v and Kolmogorov phe- 
nomenology. As long as the spectra are plotted this way and 
the scaling is correct, the curves obtained in simulations with 
different resolutions have to collapse into a single plot on 
small scales. Not only the spectrum, but any other statisti- 
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Fig. 7 Relation between the mean value of the standard deviation of density fluctuation (5p/po) and the sonic Mach number ^4K S for models 
with low resolution. Two lines show the analytical relations derived for magnetic-pressure-dominated turbulence (~ „# v ) and the gas-pressure- 



dominated turbulence (~ 



cal property of turbulence can be treated this way. When the 
resolution is increasing, the convergence is supposed to be- 
come better. This is because of the increased separation of 
scale between driving and dissipation. The optimum strat- 
egy for our MHD study is to perform the largest resolution 
simulations possible and do the scaling study. 

Figure [8] shows that the convergence study is suggestive 
of —5/3 spectral slope. The claims related to a more shallow 
slope are likely to be related to the lesser degree of locality 
of MHD turbulence, which results in the extended bottle- 



neck region ( Beresnvak and Lazarian 2010t) . The question 
definitely requires further studies, but we believe that at the 
moment there is not enough evidence to claim that GS95 
model contradicts to simulations. 



the Galactic center. We be lieve that the regime of dynamo in 
IScheko chihin et al. (2004) and the turbulence in Lazarian et al 
(2004) have similarities in terms of the density enhancement 
that are created. Although in the case of magnetic turbulence 
with suffi ciently strong mean m a gnetic field, global rever- 
sals, that iGoldreich and Sridharj (|2006) appeal to in com- 
pressing plasma, do not happen, the reversals of the mag- 
netic field direction occur in the direction perpendicular to 
the mean magnetic field. As the mean magnetic field goes 
to zero, the two regimes get indistinguishable. The regime 
of viscosity damped turbulence requires further systematic 
studies. 



3.3 Imbalanced MHD turbulence 



3.2 Extended current sheets in viscosity dominated 
turbulence 



For high Pr„, fluid Cho et al. (2002) reported a new regime 



of MHD turbulence. lLazarian et alj d2004l) showed that while 
the spectrum of volume-averaged magnetic fluctuations scales 
as Eb(k) ~ k~ l (see Figure [9), the pressure within intermit- 
tent magnetic structures increases with the decrease of the 
scale oc k and the filling factor decreases °< k . The mag- 
netic pressure compresses the gas as demonstrated in Fig- 
ure [9] More importantly, extended current sheets that natu- 
rally emerge as magnetic field fluctuates in the plane per- 
pendicular to the mean magne tic field (see Figure[9]). It was 
speculated in iLazarianl {2007) that these current sheets can 
account for the origin of the sma ll ionized and neutr a l struc- 
tures ( S INS) on AU spatial scaleslDieter et alJdl97fjl) : lHeilesl 
(ll997h : IStanimirovic et al.l (120041) . These current sheets may 
pre sent sites of cosmic ray accele ration. 



Goldreich and Sridhar (2006) appealed not to high Pr„ 



MHD turbulence per se, but to the gen eration of the mag 



netic field in the turbulent plasma (see ISchekochihin et al 



2004 to account for the high amplitude, but small scale 
fluctuations of plasma density observed in the direction of 



MHD turbulence in the presence of sources and sinks gets 
imbalanced, in the sense that the flow of energy in one di- 
rection is larger than the flow of energy in the opposite di- 
rection. Solar wind presents a vivid example of imbalanced 
turbulence with most waves near the Sun moving in the di- 
rection away from the Sun. While theories of balanced MHD 
turbulence enjoyed much attention, the theory of imbalanced 
turbulenc^H attracted less work, unfortuna tely. Important 
papers on imbalanced turbulence inclu de (see lTing et alJ 19861 
Matth aeus et alJl983 ; Biskampl2003 , and references therein). 
In terms of Solar wind observations one may mentioned 
studies bv lRoberts et alj dl987blll987al) . which showed that 
the imbalance of turbulence is not increasing contrary to 
the idealized theoretical expectations. The analy tical results 



were obtained for weak imbalanced turbulence dGaltier et al 



2002; Lit hwick and Gold reich 2003) and they are applica- 
ble in a rather narrow range of imbalance ratios. Some ear- 
lier simulations of strong imba lanced turbulence were lim 
ited to rather id ealized set ups dMaron and GoldreichfeOQl 



Cho et al. 2002), i.e. for the initial state the results of the 



simulations of strong balanced turbulence were used, but the 

13 Another name for imbalanced turbulence is a turbulence with non- 
zero cross-helicity. 
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Fig. 8 Numerical convergence of spectra in 3D reduced MHD simulations. Two upper rows are used to study convergence assuming Boldyrev 
model and two bottom rows - assuming Goldreich-Sridhar model. P(k) is so-called one-dimensional spectrum. Note that definition of dissipation 
scale r\ depends on the model, and although this difference is tiny in hyperviscous simulations Rl-5, but significa nt in viscous simu lations R6-9. 
Numerical convergence require that spectra will be similar on small scales, including the dissipation scale, see, e.g. lGotoh et al.l j2002h . As we see 
from the plots, numerical convergence is absent for Boldyrev model. For Goldreich-Sridhar model the convergence is r eached at the dissipation 
scale. Higher-resolution simulations are required to demonstrate convergence in the inertial range. Plotted with data in Beresnvak and Lazarianl 

feoioh . 




Fig. 9 Left: Filaments of density created by magnetic compression of the gas in the viscosity-damped regime of MHD turbulence. Darker regions 
correspond to higher density. The viscous damping scale l c is much larger than the current sheet thickness d. This creates large observed density 
contrasts. Center: Magnetic reversals (in the plane _L to mean (B)) that create compressions of density. Darker regions correspond to higher 
magnetic field. Right: Spectra of density and magnetic field are similar, while velocity is damped. The resistive scale in this regime is not L/Rm 
but Lffm~'/ 2 . From Beresnyak & Lazarian (preprint). 



amplitudes of waves moving in one direction were reduced. 
This did not allow making definitive conclusions about prop- 
erties of imbalanced turbulence. 



cal si mulations performed so far (Ber esnvak and Lazarian 
2010). This model can be viewed as an extension of GS95 
model into the imbalanced regime. 



Attempts to construct the mod el of stationary st rong im 



While the classic formulation of the GS95 critical bal- 



balanced turbulence were done in 



Lithwick et a 



(2008); IChandranl d2008l) ; iPerez and Boldyrev 



, (2007); BeresnvaikcBndftasazlarinn^ausalit y, is unable to describe consi stently 

d2009h. | y Lazarianl Q2008) pro- 

lowwe discuss only the model in Ber esnvak and Lazarianl posed a new way to relate parallel and perpendicular scales 
(2008) as this is the only model that agrees with numeri- of Alfvenic modes. The new condition was obtained appeal- 
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Fig. 10 Left upper: w + wavepacket, produced by cascading by w wavepacket is aligned with respect to w wavepacket, but disaligned with 
respect to the local mean field on scale X\ , by the angle 9. Left lower: the longitudinal scale A of the wavepackets, as a function of their transverse 
scale, A; Middle: the power spectrum of energies for both waves in an imbalanced forced inco mpressive 320 3 numerical simu lation. Right: the 
A (A) dependence in the same simulation, the length scales are in the units of the cube size. From lBeresnvak and La zarian (2008). 



ing to the process termed "propagation balance". In the case 
of the balanced turbulence the "old" and "new" critical bal- 
ance condition results in the same GS95 scaling. However, 
in the case of imbalanced turbulence the new formulation 
predicts smaller anisotropy for the stronger wave (which 
directly contradicts old ca usal critical balance) is consis- 
tent with simulations. The iBeresnvak and Lazarianl 020081) 
model predicts that the anisotropies of the two waves are 
different and this should be taken into account for calculat- 
ing the scattering and acceleration arising from the particle- 
wave interactions. 



4 Magnetic reconnection in astrophysical fluids 

4. 1 Problem of magnetic reconnection 

A magnetic field embedded in a perfectly conducting fluid 
is generally b elieved to preserves its topology for all time 
dParkerll 19791) . Although ionized astrophysical objects, like 
stars and galactic disks, are almost perfectly conducting, 
they show indications of changes in top ology, "magnetic re- 
connection", on dynamical tim e scales ( Parkerll970t Lovelace 



A picture of two flux tubes of different directions which 
get into contact in 3D space is a generic framework to de- 
scribe magnetic reconnection. The upper panel of Figure [TT] 
illustrates why reconnection is so slow in the textbook Sweet- 
Parker model. Indeed, the model considers magnetic fields 
that are laminar and therefore the frozen-in condition for 
magnetic field is violated only over a thin layer dominated 
by plasma resistivity. The scales over which the resistive dif- 
fusion is important are microscopic and therefore the layer 
is very thin, i.e. A <C L x , where L x is the scale at which mag- 
netic flux tubes come into contact. The latter is of the order 
of the diameter of the flux tubes and typically very large for 
astrophysical conditions. During the process of magnetic re- 
connection all the plasma and the shared magnetic flu:x0 
arriving over an astrophysical scale L x should be ejected 
through a microscopic slot of thickness A. As the ejection 
velocity of magnetized plasmas is limited by Alfven velocity 
Va, this automatically means that the velocity in the vertical 
direction, which is reconnection velocity, is much smaller 
than Va- 

We note that if magnetic reconnection is slow in some 
astrophysical environments, this automatically means that 
the results of present day numerical simulations in which 
the reconnection is inevitably fast due to numerical diffu- 



1976UPriest and Forbeskooi), Reconnection can be observed 

direct ly in the solar corona Jlnnes et all 1 9971:1 Yokovama and Shiftily do not correctly represent magnetic field dynamics in 



19951 iMasuda etalJll994 . but can also be inferred from the 



ors(Parkerlll993 


Ossendriiver 2003). Solar flares (Sturrock 


1966) and y-rav busts (Fox et al.ll2005; 


Galamaetal. 1998) 



are usually associated with magnetic reconnection. A lot of 
work in the field has concentrated on showing how recon- 
nectio n can be rapid in plasmas with very small coll i sional 
rates dShav et al.ll998llDrakel200lHDrake et al.l2006HDaughtonnj^ >n when the flare actually occurs. This bimodal charac 



these environments. If, for instance, the reconnection were 
slow in collisional media this would entail the conclusion 
that the entire crop of interstellar, protostellar and stellar 
MHD calculations would be astrophysically irrelevant. To 
make the situation more complicated, solar flares demon- 
strate both periods of slow reconnection when the accumu- 
lation of magnetic flux takes place and periods of fast recon- 



2006), which substantially constrains astrophysical applica- 



tions of the corresponding reconnection models. The dis- 
regard of pre-existing turbulence is another shortcoming of 
these models. 



14 FigurefTJIpresents a cross section of the 3D reconnection layer. A 
guide field is present in the generic 3D configurations of reconnecting 
magnetic flux tubes. 
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ter is another property that a successful model of reconnec- 
tion should address. 

In this section we focus our attention on the reconnection 
that takes place in turbulent astrophysical flui ds. The recon- 



and particle acceleration throughout the univer se (Ide Gouveia Dal Pino and 
20031 120051 lLazarianll2005l brake et alJl2006h . 



4.2 R ate of reconnection 
nection in such environments was described in Lazarian and Vishniac 



dl999h . The latter model identified magnetic field wander- 
ing as the key process that induces fast, i.e. independent of 
fluid resistivity, magnetic reconnection. The predicted de- 
pendences of magnetic reconnection on the properties of 
turbulence, i.e. on the intensity and inje ction scale of turbu- 
lence , have been successfully tested in iKowal et al.1 (12009 , 
2012b . The deep consequence of the model is the violation 
of the magnetic flux freezing in turbulent fluids (e .g. the ex- 
plicit claim of this in lVishniac and Laza rian 1999). The con- 
nection of the LV99 model with more recent ideas on the 
properti es of magnetic fiel d in turbulent plasmas are ana- 
lyzed in lEvink et al.l (1201 II) . 



While the discussion of possible ways how turbulence 
can enhance reconnection rates was not unprecedented, the 
LV99 model was radically different from what was proposed 
earlier. For in stance, the closest to th e spirit of LV 99, are 
two papers by Mat thaeus and Lamkinl d 19851 Il986h . where 
the authors performed 2D numerical simulations of turbu- 
lence and provided arguments in favor of magnetic recon- 
nection getting fast. However, the physics of the processes 
that they considered was very different from that in LV99. 
For instance, the key process of field wandering of the LV99 
model has n ot been considered by iMatthaeus and Lamkin 



( 1985 , 1986b . On the contrary, the components of their ap- 



proach, e.g. X-point and possible effects of heating and com- 
pressibilit y are no t ingredient s of the LV99 model . Other pa- 
pers, e.g. ISpeiserjdl970l) and lJacobson an d Moses ( 1984J) are 
even more distant in terms of the effects that they explored, 
namely, they studied the changes of the microscopic proper- 
ties of the plasma induced by turbulence and considered how 
these changes can accelerate magnetic reconnection. At the 
same time, LV99 shows that the microscopic plasma prop- 
erties are irrelevant for their model of reconnection. 

Below (Sj5]l we present numerical evidence, based on three 
dimensional simulations, that reconnection in a turbulent fluid 
occurs at a speed comparable to the rms velocity of the tur- 
bulence, regardless of either the value of the resistivity or 
degree of collisionality. In particular, this is true for turbu- 
lent pressures much weaker than the magnetic field pressure 
so that the magnetic field lines are only slightly bent by the 
turbulence. These results are consistent with the proposal by 
LV99 that reconnection is controlled by the stochastic diffu- 
sion of magnetic field lines, which produces a broad outflow 
of plasma from the reconnection zone. This implies that re- 
connection in a turbulent fluid typically takes place in ap- 
proximately a single eddy tu rnover time, with broad impli- 
cations for dynamo activity dParkeill 1 97d 1 1 9931: Istixll2000h 



Astrophysic al plasmas are often highly ionized and highly 
magnetized ( Parker 1970l) . The evolution of the magnetic 
field in a highly conducting fluid can be described by a sim- 
ple version of the induction equation 



= Vx(vxB-t]VxB) 



(7) 



~di 

where B is the magnetic field, v is the velocity field, and rj 
is the resistivity coefficient. Under most circumstances this 
is adequate for discussing the evolution of magnetic field 
in an astrophysical plasma. When the dissipative term on 
the right hand side is small, as is implied by simple dimen- 
sional estimates, the magnetic flux through any fluid ele- 
ment is constant in time and the field topology is an in- 
variant of motion. On the other hand, reconnection is ob- 



served in the solar corona and chromosphere (llnnes etal, 

ll997llYokovama and Shibatall995llMasuda et aljl994 ; Ciaravella and Rav 



2008), its presence is re quired to explain d ynamo action in 
stars and galactic disks ( Parker|[l970l 1993 ). and the violent 



relaxation of magnetic fields following a change in topology 
is a prime candidate for the acceleration of high energy parti- 
cl es (de Gouveia Dal Pino and Lazarianl 2003. 2005; Lazarian 



2005 



Drake et al.ll2006l lLazarianll2009l: brake et alE oiO) 



in the universe. Quantitative general estimates for the speed 
of reconnection start with two adjacent volumes with differ- 
ent la rge scale magnetic field directions ( Sweet! 1958 ; Parker 



1957). 



The speed of reconnection, i.e. the speed at which in- 
flowing magnetic field is annihilated by Ohmic dissipation, 
is roughly tj /A, where A is the width of the transition zone 
(see Figure [TTJ. Since the entrained plasma follows the lo- 
cal field lines, and exits through the edges of the current 
sheet at roughly the Alfven speed, Va, the resulting recon- 
nection speed is a tiny fraction of the Alfven speed, Va = 
B / (4-np) 1 ! 2 where L is the length of the current sheet. When 
the current sheet is long and the reconnection speed is slow 
this is referred to as Sweet-Parker reconnection. Observa- 
tions require a speed close to Va, so this expression implies 
that L ~ A, i.e. that the magnetic field lines reconnect in an 
"X point". 

The f i rst mo del with a stable X point was proposed by 
Petschekl (1964). In this case the reconnection speed may 



have little or no dependence on the resistivity. The X point 
configuration is known to be unstable to collapse into a sheet 
in the MHD regime (see Biskampl 1996 ). but in a collision- 
less plasma it can be maintaine d through coupling to a dis- 
persive plasma mode ( Sturrocklll966h . This leads to fast re- 
connection, but with important limitations. This process has 
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a limited astrophysical applicability as it cannot be impor- 



tant fo r most phases of the interstellar medium (see Dra ine and 
19981 for a list of the idealized phases), not to speak about 
dense plasmas, such as stellar interiors and the denser parts 
of accretion disks. In addition, it can only work if the mag- 
netic fields are not wound around each other, producing a 
saddle shaped current sheet. In that case the energy required 
to open up an X point is prohibitive. The saddle current sheet 
is generic for not parallel flux tubes trying to pass through 
each other. If such a passage is seriously constrained, the 
magnetized highly conducting astrophysical fluids should 
behave more like Jello rather than normal fluids^ 

Finally, the traditional reconnection setup does not in- 
clude ubiquitous astrophysical turbulenc J^lfsee Armstrong et 



Lazarian 




Sweet-Parker model 




1995llE lmegree n and ScalobOQ^Iiv lcKee an d Ostrikerl2007l: 
Lazar ian 2009; Chepurnov and Lazarian 2010). Fortunately, 
this approach provides another way of accelerating recon- 
nection. Indeed, an alternative approach is to consider ways 
to decouple the width of the plasma outflow region from 
A . The plasma is constrained to move along magnetic field 
lines, but not necessarily in the direction of the mean mag- 
netic field. In a turbulent medium the two are decoupled, 
and fluid elements that have some small initial separation 
will be separated by a large eddy scale or more after moving 
the length of the current sheet. As long as this separation is 
larger than the width of the current sheet, the result will not 
depend on T] . 

LV99 introduced a model that included the effects of 
magnetic field line wandering (see Figure [TTJ. The model 
relies on the nature of three-dimensional magnetic field wan- 
dering in turbulence. The effects of compressibility and heat- 
ing, which were thought to be impor tant in the earlier studies 
(see, for ex ample, the discussion in iMatthaeus and Lamkin 
1985i Il986h . do not play the role for the LV99 model ei- 
ther. The model is applicable to any weakly perturbed mag- 
netized fluid, irrespectively, o f the degree of p lasma being 
collisional or collisionless (cf. Sha v et al ,1ll998h . 

Two effects are the most important for understanding of 
the nature of reconnection in LV99. First of all, in three 
dimensions bundles of magnetic field lines can enter the 
reconnection region and reconnection there independently 
(see Figure [TTJ, which is in contrast to two dimensional 
picture where in Sweet-Parker reconnection the process is 
artificially constrained. Then, the nature of magnetic field 
stochasticity and therefore magnetic field wandering (which 
determines the outflow thickness, as illustrated in Figure fTTt 

The idea of felt or Jello-type structure of magnetic fields in inter- 
stellar medium was advocated by Don Cox (private communication). 

16 The set ups where in stabilities play important role include 
IShimizu et alj {2009b 2009a). For sufficiently large resolution of sim- 
ulations those set-ups are expected to demonstrate turbulence. Turbu- 
lence initiation is also e xpected in the presence of plasmoid ejection 
(Shibata and Tanuma 2001). Numerical viscosity constrains our ability 
to sustain turbulence via reconnection, however. 




Turbulent model 



blow up 

Fig. 11 Upper plot: Sweet-Parker model of reconnection. The outflow 
is limited by a thin slot A, which is determined by Ohmic diffusivity. 
The other scale is an astrophysical scale L>4. Middle plot: Recon- 
nection of weakly stochastic magnetic field according to LV99. The 
model that accounts for the stochasticity of magnetic field lines. The 
outflow is limited by the diffusion of magnetic field lines, which de- 
pends on field line stochasticity. Low plot: An individual small scale 
reconnection region. The reconnection over small patches of magnetic 
field determines the local reconnection rate. The global reconnection 
rate i s substantially larger a s many independent patches come together. 
From[Lazarian et al. (2 0041) . 



is very different in 2D and the real 3D world (LV99). In other 
words, removing artificial constraints on the dimensionality 
of the reconnection region and the magnetic field being ab- 
solutely straight, LV99 explores the real-world astrophysical 
reconnection. 

The scaling relations for Alfvenic turbulence allow us to 
calculate the rate of magnetic field spreading. This is the ef- 
fect known to be important for the cosmic ray propagation 
Jokipiil (1 197 3b and LV99 quantified it for the case of GS95 



turbulencd 17 l A bundle of field lines confined within a region 
of width y at some particular point will spread out perpen- 
dicular to the mean magnetic field direction as one moves in 
either direction following the local magnetic field lines. The 
rate of field line diffusion is given approximately by 



dx 



A,, 



(8) 



This treatment was used later for deali ng with different prob- 
lems, f rom heat propagation i n pla smas iNaravan and Medvedevl 
(2001); lLazarian and Pogosvan (2006) to cosmic ray propagation 



lYan and Laz arian (2008). 



18 



where Am 1 ~ l» , t\\ is the parallel scale chosen so that the 

corresponding vertical scale, £j_, is ~ (y 2 ) 1 ^ 2 , and x is the 
distance along an axis parallel to the mean magnetic field. 
Therefore, using equation (0 one gets 



tlx 



where we have substituted (y 2 ) 1 / 2 for £j_. This expression 
for the diffusion coefficient will only apply when y is small 
enough for us to use the strong turbulence scaling relations, 
or in other words when (y 2 ) < L 2 (ul/va) 4 - Larger bundles 
will diffuse at a maximum rate of L;(ml/va) 4 - For (y 2 ) small 
equation (0 implies that a given field line will wander per- 
pendicular to the mean field line direction by an average 
amount 

T 3/2 /.. x2 



(/) 1/2 -^T7l(- 



(10) 



in a distance x. The fact that the rms perpendicular displace- 
ment grows faster than x is significant. It implies that if we 
consider a reconnection zone, a given magnetic flux element 
that wanders out of the zone has only a small probability of 
wandering back into it. 

The advantage of the classical Sweet-Parker scheme of 
reconnection is that it naturally follows from the idea of 
Ohmic diffusion. Indeed, mass conservation requires that the 
inflow of matter through the scale of the contact region L x 
be equal to the outflow of matter through the diffusion layer 
A, i.e. 

A 

Vrec = VA — - (11) 

At 

The mean-square vertical distance that a magnetic field- 
line can diffuse by resistivity in time t is 



(y 2 (t))~Xt. 



(12) 



The field lines are advected out of the sides of the recon- 
nection layer of length L x at a velocity of order va . Thus, 
the time that the lines can spend in the resistive layer is the 
Alfven crossing time t& = Lx/va- Thus, field lines can only 
be merged that are separated by a distance 

4 = yj (y 2 (t A )) ~ V^a~ = Lx/Vs, (13) 

where S is Lundquist number. Combining Eqs. (fTTT l and ( fT3l 
one gets the famous Sweet-Parker reconnection rate, v rec — 
va/VS. 

In LV99 magnetic field wandering determines the scale 
of the outflow A (see FigurefTTI). Using expressions from the 
earlier section one can obtain (LV99): 

1/2 / t \ 1/2] 



Vrec < va min 



{ul/VaY 



(14) 



This limit on the reconnection speed is fast, both in the sense 
that it does not depend on the resistivity, and in the sense that 



it represents a large fraction of the Alfven speed. To prove 
that Eq. (TBI indeed constitutes the reconnection rate LV99 
goes through a thorough job of considering all other possible 
bottlenecks for the reconnection and shows that they provide 
higher reconnection speed. 

Below we provide a new derivation of the LV99 recon- 
nection rates which makes apparent that the LV99 model is 
a natural generalization of the laminar Sweet-Parker model 
to flows with back ground turbulence. The new argument in 
Evink et alJ (1201 if) is based on the concept of Richardson 
diffusion. It is known in hydrodynamic turbulence that the 
combination of small scale diffusion and large scale shear 
gives rise to Richardson diffusion, where the mean square 
separation between two particles grows as f 3 once the rms 
separation exceeds the viscous damping scale. A similar phe- 
nomenon occurs in MHD turbulence. In both cases the sep- 
aration at late times is independent of the microscopic trans- 
port coefficients. Although the plasma is constrained to move 
along magnetic field lines, the combination of turbulence 
and ohmic dissipation produces a macroscopic region of points 
that are "downstream" from the same initial volume, even in 
the limit of vanishing resistivity. 

Richardson diffusion (see iKupiainenlEobH) implies the 



mean squared separation of particles (\x\ (f) — x%(f )| 2 ) rs £f 3 , 
where t is time, e is the energy cascading rate and (•) de- 
notes an ensemble averaging. For subAlfvenic turbulence 
e u A L / (v^L,-) (see LV99) and therefore analogously to Eq. 
( fTST l one can write 



! L(L/L i ) 1 / 2 (« i /V A ) 2 



(15) 



where it is assumed that L < L,. Combining Eqs. ( TTTb and 
( TT3T > one gets 

V r ec, L V99 « VA(L/Li)^ 2 (u L /V A ) 2 . (16) 

in the limit of L < L,-. Analogous considerations allow to 
recover the LV99 expression for L > Li, which differs from 
Eq. (IToT l by the change of the power 1/2 to —1/2. These 
results coincide with those given by Eq. (fT4l . 

It is important to stress that Richardson diffusion ulti- 
mately leads to diffusion over the entire width of large scale 
eddies once the plasma has moved the length of one such 
eddy. The precise scaling exponents for the turbulent cas- 
cade does not affect this result, and all of the alternative 
scalings considered in LV99 yield the same behavior. 

Other forms of reconnection rate are also useful. For in- 
stance, one can take into account that the velocity of the tran- 
sition to the strong turbulent regime is u tur i, strong ~ VaM\. 
Then 



V re 



M-turb,strong (} / L, 



A 1/2 



(17) 



shows that the reconnection rate scales with the velocity 
of the strong turbulence cascade. Similarly, it is useful to 
rewrite this expression in terms of the power injection rate 
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Pi„j. As the perturbations on the injection scale of turbu- 
lence are assumed to have velocities U[ < Va, the turbulence 
is weak at large scales. Therefore, the relation between the 
power and the injection velocities are different from the usual 
Kolmogorov estimate, namely, in the case of the weak tur- 
bulence P ~ u\ /(IV A ) (LV99). Thus we get, 



Vr, 



P A 1/2 



(18) 



where I is the length of the turbulent eddies parallel to the 
large scale magnetic field lines as well as the injection scale. 

If turbulence is superAlfvenic, LV99 predicts that the 
diffusion of magnetic field happens with the turbulent ve- 
locity. This is possible as the subAlfvenic reconnection that 
happens below the scale I a given by Eq. © does not present 
the bottleneck for the reconnection at the larger scales. 

The LV99 reconnection velocity (e.g. Eq. ( TT8l ) does not 
depend on resistivity or plasma effects. Therefore for suffi- 
ciently high level of turbulence we expect both collisionless 
and collisional fluids to reconnect at the same rate. 



4.3 Model of flares 

If turbulence can drive reconnection, which in turn trans- 
forms magnetic energy into kinetic energy, then it seems 
appropriate to wonder if the process can be self-sustaining. 
That is, given a very small level of ambient turbulence, how 
likely is it that reconnection will speed up as it progresses, 
without any further input from the surrounding medium? 
The corresponding problem was in LV99. Below we follow 
simple arguments provided in lLazarian and Vishniac ( 2009 ), 
which clarify the related physics by presenting an idealized 
model of a reconnection flare. 

Let' s consider a reconnection region of length L and thick- 
ness A . The thickness is determined by the diffusion of field 
lines, which is in turn determined by the strength of the tur- 
bulence in the volume. Reconnection will allow the mag- 
netic field to relax, creating a bulk flow. However, since 
stochastic reconnection is expected to proceed unevenly, with 
large variations in the current sheet, we can expect that some 
unknown fraction of this energy will be deposited inhomo- 
geneously, generating waves and adding energy to the local 
turbulent cascade. We take the plasma density to be approx- 
imately uniform so that the Alfven speed and the magnetic 
field strength are interchangeable. The nonlinear dissipation 
rate for waves is 



k 2 v 2 

K ± v wave £ 



hV A 



> K i_ v turb^.L,turb 



(19) 



where the first rate is the self-interaction rate for the waves 
and the seco nd is the dissipation rate by t he ambient tur- 
bulence (see iBeresnvak and Lazari an 2008). The important 



point here is that k± for the waves falls somewhere in the in- 
ertial range of the strong turbulence. Eddies at that wavenum- 
ber will disrupt the waves in one eddy turnover time, which 
is necessarily less than L/Va- The bulk of the wave energy 
will go into the tubulent cascade before escaping from the 
reconnection zone. (This zone will radiate waves, for the 
same reason that turbulence in general radiates waves, but 
it will not significantly impact that energy budget of the re- 
connection region.) 

We can therefore simplify our model for the energy bud- 
get in the reconnection zone by assuming that some fraction 
£ of the energy liberated by stochastic reconnection is fed 
into the local turbulent cascade. The evolution of the turbu- 
lent energy density per area is 



(Av 



^ V — ' turb J ~-a i k t- • niro — 

where the loss term covers both the local dissipation of tur- 
bulent energy, and its advection out of the reconnection zone. 
Since V rec ~ v tur j, and A ~ L(y tur blV A ), we can rewrite this 
by defining J{, A = v tur b /Va and T = L/Va so that 

d 



= SVjtVre, 



V A 



l turb L 



(20) 



dx 

If e is a constant then 



■■V A e 



1/2 



1-1- 



,-2t/3 



-,1/2 



(21) 



(22) 



This implies that the reconnection rate rises to £ 1 I 2 Va is a 
time comparable to the ejection time from the reconnection 
region (~ L/Va)- Given that reconnection events in the solar 
corona seem to be episodic, with longer periods of quies- 
cence, this implies that either e is very small, for example - 
dependent on the ratio of the thickness of the current sheet 
to A , or is a steep function of ^a ■ If it scales as JK A to some 
power greater than two then initial conditions dominate the 
early time evolution. 

An alternative route by which stochastic reconnection 
might be self-sustaining would be in the context of a series 
of topological knots in the magnetic field, each of which 
is undergoing reconnection. Now the problem is sensitive 
to geometry. Let's assume that as each knot undergoes re- 
connection it releases a characteristic energy into a volume 
which has the same linear dimension as the distance to the 
next knot. The density of the energy input into this volume 
is roughly e V^v tur i, /L, where e is the efficiency with which 
the magnetic energy is transformed into turbulent energy. 
We have 

,,'3 



eVi 



2 v turb v ' 



(23) 



L L k : 

where L k is the distance between knots and v' is the turbulent 
velocity created by the reconnection of the first knot. This 
process will proceed explosively if v' > v, u ,-i or 

V 2 L k e > v 2 urb L. (24) 
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This condition is almost trivial to fulfill. The bulk motions 
created by reconnection will unavoidably generate signifi- 
cant turbulence as they interact with their surrounding, so e 
should be of order unity. Moreover the length of any cur- 
rent sheet should be at most comparable to the distance to 
the nearest distinct magnetic knot. The implication is that 
each magnetic reconnection event will set off its neighbors, 
boosting their reconnection rates from v tU rb, set by the envi- 
ronment, to e l l 2 VA{Lk/L) l l 2 (as long as this is less than Va). 
The process will take a time comparable to L/v tur b to begin, 
but once initiated will propagate through the medium with 
a speed comparable to speed of reconnection in the individ- 
ual knots. In a more realistic situation, the net effect will be 
a kind of modified sandpile model for magnetic reconnec- 
tion in the solar corona and chromosphere. As the density of 
knots increases, and the energy available through magnetic 
reconnection increases, the chance of a successfully propa- 
gating reconnection front will increase. 

The flare of reconnection presents a reconnection insta- 
bility. In the absence of external turbulence the original out- 
flow, e.g. originated through tearing instability (see Loureiro 
et al. 2007, Bhattacharjee et al. 2009) gets turbulent and trig- 
gers the mechanism described above. This shows that the 
tearing and turbulent mechanisms may be complementary^ 



5 Testing of LV99 model 

Here we describe the results of a series of three dimensional 
numerical simulations aimed at adding turbulence to the sim- 
plest reconnection scenario and testing equation (f 1 St . We di- 
vide the domain in the middle into two equal regions along 
the vertical direction (along the Y axis in our case) with 
equal magnetic field strength, but opposite sign of its X com- 
ponent. In addition, a constant guide field is added in the Z 
direction. The domain is periodic along the guide field (i.e. 
the Z axis) and open in the other directions (the X and Y 
axis). The external gas pressure is uniform and the magnetic 
fields at the top and bottom of the box are taken to be the 
specified external fields plus small perturbations to allow for 
outgoing waves. The grid size in the simulations varied from 
256x512x256 to 512x1028x512 so that the top and bottom 
of the box are far away from the current sheet and the region 
of driven turbulence around it. At the sides of the box where 
outflow is expected the derivatives of the dynamical vari- 
ables are set to zero. A complete d escription of the nu mer- 
ical methodology can be found in iKowal et al.1 d2009h . All 
simulations are allowed to evolve for seven Alfven cross- 
ing times without turbulent forcing. During this time they 

18 When turbulence develops the LV99 mechanism can provide much 
faster reconnection compared to tearing and tearing may become a sub- 
dominant process. In fact, emerging turbulence may suppress the tear- 
ing instability. 



develop the expected Sweet-Parker current sheet configura- 
tion with slow reconnection. Subsequently the isotropic tur- 
bulent forcing is turned on inside a volume centered in the 
midplane (in the XZ plane) of the simulation box and ex- 
tending outwards by a quarter of the box size. The turbu- 
lence reaches its full amplitude around eight crossing times 
and is stationary thereafter. 

The speed of reconnection in three dimensions can be 
hard to define without explicit evaluation of the magnetic 
field topology. However, in this simple case we can define it 
as the rate at which the x component of the magnetic field 
disappears. More precisely, we consider a yz slice of the 
simulation, passing through the center of the box. The rate 
of change of the area integral of — B x — is its flux across 
the boundaries of the box minus the rate at which flux is 
annihilated throug h reconnection (see more discussion in 
Kowal et alj|2009l) 



dt(j \B x \dzdyj = j sign(B x )Edl - 2V rec B x , ext L z (25) 

where electric field is E = v x B — rjj, B x ext is the absolute 
value of B x far from the current sheet and L z is the width 
of the box in the z direction. This follows from the induc- 
tion equation under the assumption that the turbulence is 
weak to lead to local field reversals and that the stresses at 
the boundaries are weak to produce significant field bending 
there. In other words, fields in the x direction are advected 
through the top and bottom of the box, and disappear only 
through reconnection. Since periodic boundary conditions 
are adopted in the z direction the boundary integral on the 
right hand side is only taken over the top and bottom of the 
box. By design this definition includes contributions to the 
reconnection speed from contracting loops, where Ohmic 
reconnection has occurred elsewhere in the box and \B X \ de- 
creases as the end of a reconnected loop is pulled through 
the plane of integration. It is worth noting that this estimate 
is roughly consistent with simply measuring the average in- 
flux of magnetic field lines through the top and bottom of 
the computational box and equating the mean inflow veloc- 
ity with the reconnection speed. Following equation d25l l we 
can evaluate the reconnection speed for varying strengths 
and scales of turbulence and varying resistivity. 

In Figure[13]we see the results for varying amounts of in- 
put power, for fixed resistivity and injection scale as well as 
for the case of no turbulence at all. The line drawn through 
the simulation points is for the predicted scaling with the 
square root of the input power. The agreement between equa- 
tion dT~8b and Figure[l3]is encouraging but does not address 
the most important aspect of stochastic reconnection, i.e. its 
insensitivity to f] . 

In Figure[l4]we plot the results for fixed input power and 
scale, while varying the background resistivity. In this case 
rj is taken to be uniform, except near the edges of the com- 
putational grid where it falls to zero over five grid points. 
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Fig. 12 Visualization of reconnection simulations in Kow al et"aTI 120121) . Left panel: Magnetic field in the reconnection region. Large perturbations 
of magnetic field lines arise from reconnection rather than driving; the latter is subAlfvenic. The color corresponds to the direction of magnetic 
lines. Central panel: Current intensity and magnetic field configuration during stochastic reconnection. We show a slice through the middle of the 
computational box in the xy plane after twelve dynamical times for a typical run. The guide field is perpendicular to the page. The intensity and 
direction of the magnetic field is represented by the length and direction of the arrows. The color bar gives the intensity of the current. The reversal 
in B x is confined to the vicinity of y=0 but the current sheet is strongly disordered with features that extend far from the zone of reversal. Right 
panel: Representation of the magnetic field in the reconnection zone with textures. 



This was done to eliminate edge effects for large values of 
the resistivity. We see from the Figure [14] that the points for 
laminar reconnection scale as ^/rf, the expected scaling for 
Sweet-Parker reconnection. In contrast, the points for recon- 
nection in a turbulent medium do not depend on the resistiv- 
ity at all. In summary, we have tested the model of stochas- 
tic reconnection in a simple geometry meant to approximate 
the circumstances of generic magnetic reconnection in the 
universe. Our results are consistent with the mechanism de- 
scribed by LV99. The implication is that turbulent fluids in 
the universe including the interstellar medium, the convec- 
tion zones of stars, and accretion disks, have reconnection 
speeds close to the local turbulent velocity, regardless of the 
local value of resistivity. Magnetic fields in turbulent fluids 
can change their topology on a dynamical time scale. 

Figure[T5]shows results of numerical experiments in which 
the dependence of the reconnection on the anomalous resis- 
tivity was studied. The anomalous resistivity increases ef- 
fective resistivity for high current densities. It is frequently 
used as a proxy for plasma effects, e.g. collisionless effects 
in reconnection. While this type of resitistivity enhances the 
local speed of indi vidual reconnection events, the study in 
Kowal et al.1 (120091) testifies that the total reconnection rate 
does not change. 

Any numerical study has to address the issue of the pos- 
sible numerical effects on the results. We show the depen- 
dence of the reconnection rate on the numerical resolution in 
Figure[l6] The reconnection rate increases with the increase 
of the resolution, which testifies that the fast reconnection 
is not due to numerical effects. Indeed, higher numerical re- 
connection is expected for lower resolution simulations. 



Finally, it is important to give a few words in relation to 
our turbulence driving. We drive our turbulence solenoidally 
to minimize the effects of compression, which does not play 
a role in LV99 model. The turbulence driven in the volume 
around the reconnection layer corresponds to the case of as- 
trophysical turbulence, which is also volume-driven. On the 
contrary, the case of the turbulence driven at the box bound- 
aries would produce spatially inhomogeneous imbalanced 
turbulence for which we do not ha ve analytical predictions 
(see d iscussion of such turbulence in lBeresnvak and Lazarian 
2009). We stress, that it is not the shear size of our numerical 
simulations, but the correspondence of the observed scalings 
to those predicted in LV99 that allows us to claim that we 
proved that the 3D reconnection is fast in the presence of 
turbulence. 



6 Turbulence and cosmic rays acceleration 

6.1 Shock acceleration in the presence of turbulence 

Cosmic rays (CRs), relativistic charged particles with ener- 
gies 10 8 — 10 22 e V, constitute an es sential part of astrophys- 
ical systems (see lSchlickeiserll2003l) . The origin of CRs has 
been a su bject of debate from the begin ning of research in 
the field ( Ginzburg and Svrovatskiil["l964 ). By now, it is ac- 
cepted that galactic CRs at least up to the "knee" in the spec- 
trum, just above 10 15 eV, are most likely generated primarily 
by strong supernova shocks. However, it is easy to show that 
typical interstellar upstream magnetic field of around 5/iG 
is too weak to provide the return of the cosmic rays with 
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Fig. 13 Reconnection speed versus input power Pi n j for the driven 
turbulence. Upper plot: Variations of reconnection speed in time for 
different levels of turbulent driving. Lower Plot: Reconnection speed, 
plotted against the input power for an injection wavenumber equal to 8 
(i.e. a wavelength equal to one eighth of the box size) and a resistivity 
v„. The dashed line is a fit to the predicted dependence of P 1 / 2 . The 
horizontal line shows the laminar reconnection rates for each of the 
simulations before the turbulent forcing started. Here the uncertainty 
in the time averages are indicated by the si ze of the symbols a nd the 
variances are shown by the error bars. From lKowal et al. (2009). 
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Fig. 14 Reconnection speed versus resistivity. We show the reconnec- 
tion speed plotted against the uniform resistivity of the simulation for 
an injection wavenumber of 8 and an injected power of one. We in- 
clude both the laminar reconnection speeds, using the hollow symbols, 
fit to the expected dependence of r/„, and the stochastic reconnection 
speeds, using the filled symbols. As before the symbol sizes indicate 
the uncertainty in the average reconnection speeds and the error bars 
indicate the variance. B x = 1 and simulations with large, B 7 = 1 , and 
small, B z = 0.1, guide fields are shown. FromK<walj;tal. (2009). 
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energies ~ 10 15 GeV back to the shock to continue the ac- 
celeration. 

To overcome the problem it is important to increase the 
value of the magnetic field both in the preshock and post- 
shock regions. The latter problem has a straightforward so- 
lution, as postshock gas is known to be turbulent. Turbu- 



Fig. 15 Effect of anomalous r esistivity on the rec onnection speed. No 
dependence is observed. From lKowal et al. (2009). 



2007tlZirakashvili et alJ2008[|Riquelme and Spifkovskvl2009l) 



with many researchers pointing out that the instability may 
become much suppressed. In addition, even at its linear regime 
the instability is lik ely to be suppressed in the presence of 



lence can amplify the magnetic field ( see Kazantsev 1968t the instability is lik ely to be suppressed in the presence ot 
Kulsrud and Andersonll992Hcho et alj2009l:lBeresnvak et alj ambient turbulence jYan and Lazarianl2002l:lFar mer and Goldreich 

2004t iBeresnvak and Lazarianl 120081) . In th i s situa tion, the 



201 1) and therefore an enhan ced magnetic field in the post - 



current-driven instability proposed by iBeil J2004 became 
the central idea being discussed in order to solve the prob- 
lem of the acceleration of high energy cosmic rays. 

Beresnyak, Jones & Lazarian (2009, henceforth BJL09) 
suggested an alternative model for enhancing the magnetic 
field in the preshock region. This model is based on the 
mag netic field generation within the shock precursor. The 
Blasi and Amatol2008tlLucek and BeHl2000tlDiamond and Malknreb ursor originates from cosmic rays that can be returned 



shock region is expected (see iGiacalone and Jok ipii 2007). 
A more challenging issue is to have the magnetic field am- 
plified in preshock plasmas. One known way is to appeal to 
the streaming instability. The latter evolves to the nonlinear 
stage with 8B ^ Bo making the original classical treatment 
of the instability not applicable. What happens in the non- 
linear regime has been discussed a lot in recent years (e.g. 
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Fig. 16 Dependence of the reconnection rate on the numerical resolu- 
tion. If the fast reconnection were due to yet unclear numerical effects 
on small scales, we would expect to see the increase of the reconnec- 
tion rate with the decrease of the numerical box. If anything, the actual 
dependence of the reconnection rate on the box size shows the opposite 
dependence. From Kowal et al. ( 2 0121) . 
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Fig. 17 Magnetic field spectrum (dashed), generated by a small-scale 
dynamo induced by solenoidal velocity motions (solid). L* is an 
equipartition scale of magnetic and kinetic motions, it plays a central 
role in particle sc attering. Upper p anel: magnetic and velocity fields 
from simulations ( ICho et alj |2009) with different dashed lines corre- 
sponding to magnetic spectra at different times. From BIL09. 



back to the shock by the original interstellar magnetic field. 
The fluid is stirred within the precursor by the combina- 
tion of fluid density inhomogeneities and CR pressure. This 
drives vorticity and the precursor turbulence. The latter am- 
plifies magnetic energy through a turbulent dynamo. In other 
words, B JL09 claims that both the enhancement of magnetic 
field in the preshock and postshock regions are due to turbu- 
lence draining energy from the shock. 

The upstream diffusion of CRs leading to precursor for- 
mation is associated with a n umber of possibilities, such as 
the Drury acoustic ins tability dDrurvlll984 ; Dorfi and Drurv 
19851: lKanget~aTll992h . which is the enhancement of com- 
pressible perturbations by the CR pressure gradient. Such in- 
stabilities, however, will only operate on perturbations dur- 
ing the limited time that a fluid element crosses the precur- 
sor. This crossing time is rather short T £ « l p /U s where l p is 
the precursor thickness and U s is the shock velocity and the 
growth of the density perturbation is constrained. 

The difference of the velocities of the fluid elements 
crossing the shock arises from the difference of the iner- 
tia of dense and rarefied gas entering the shock. When an 
inhomogeneous fluid enters the precursor with speed uq it 
gets decelerated by the CR pressure gradient until it reaches 
the speed u\ at the dissipative shock front. Let us denote 
A s (uq — mi) as the difference between the ballistic velocity 
of the high-density region and the full decelerated one of 
the low density regions. The existence of such regions in in- 
terstellar medium or solar wind plasma, is expected to arise 
from the pre-existing turbulence not related to the shock (see 
SjlJ, Figure [TJ. The interaction of the precursor with the den- 
sity inhomogeneities is expected to result in intensive turbu- 
lence in the precursor. The turbulence in conducting fluid is 
known to amplify magnetic fields as depicted in Figure [T71 



After extremely short kinematic stage, the energy is grow- 
ing linearly with time dEg/dt s» a^E, where e ~ pu%/L is 
the e nergy transfer r ate, and a c / is the efficiency factor ~ 
0.06 dCho et al.ll 2009). The turbulent dynamo gets magnetic 
energy in the equipartition with the kinetic energy at a scale 
L*. This scale grows with time. On scales smaller than L* 
MHD turbulence exists. On scales larger than L* , the turbu- 
lence is superAlfvenic and velocity has a Kolmogorov spec- 
trum. 

As the crossing time of the precursor x c is short the gen- 
eration of magnetic field in the precursors does not get into 
the saturation regime. The system of equations suggested in 
BJL09 for the evolutions of magnetic field and L* is as fol- 
lows: 



8B 2 {L*, Xl ) = inA d ex(x 1 ); 

SB* f L*(xi) \ 1/3 

~ Us \ L 

. . [* dx 
t(*i) = / "Tv 

L*(x 1 ) = (2A rf «. ( T(x 1 )) 3 / 2 L- 1 /2. 



(26) 
(27) 

(28) 
(29) 



This approach was used to calculate the diffusion coeffi- 
cients shown in Figure ITU 



It was shown by Bell (2004) that magnetic fields in the 
preshock region could be amplified significantly compared 
to the background field. The model is based on a current- 
driven instability. It was successfully tested numerically in a 
numberof papers ( Vladimirov et al.l2006l: lzirakashvili et al.1 
20081 iRiquelme and Spitkovskvll2009l) . In iBelll J2004 the 
growth of magnetic field at the linear stage is fastest on the 
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instability arises only from those CRs tha t are able to freely 
stream. The efficiency can be reduced and Belli ( 2004) insta- 
bility wins if A s parameter is much less than unity. Further 
research should clarify the relative importance of the two 
ways of generating magnetic fields in the precursor. BJL09 
estimate that the maximum energy that can be obtained by 
the acceleration of the shocks with the small-scale dynamo 
in the precursor is ~ 10 16 eV, which is consistent with the 
galactic cosmic ray measurements. 



6.2 Second order Fermi acceleration and turbulence 

To deal with the problem of propagation and acceleration 
of CRs the so-called diffusive approximation is frequently 
used. Within this approximation it is assumed that the par- 
ticle scatter or gain energy in small steps and the dynamics 
is averaged to obtain the spatial diffusion coefficient, Dxx, 
and the momentum diffusion coefficient, D pp . The result- 
ing advection-diffusion equation for the evolution of quasi- 
isotropic CR distribution function / is 



Fig. 18 Diffusion coefficients. (1) is low-energy scattering, which de- 
pends on the properties of small-scale MHD turbulence (dashed: fast 
modes are present; dotted: fast modes are absent); (2) is medium- 
energy scattering, by the magnetic fields generated by turbulent dy- 
namo; (3) is high-energy scattering. From BJL09. 



characteristic scale that is determined by the initial field Bq, 
and the current 

cB Q 



I = l/kc 



(30) 



where j ( [ arises from high-energy CRs that correspond to 
qB$/ pc << k c . The linear growth rate depends only on the 
current according to the relation 



Jd 



(31) 



BLJ09 demonstrated that, conservatively assuming that 
the linear growth rate is applicable to the current-driven in- 
stability is in the stage with 8B ~ Bo, one still gets that the 
dynamo generation of magnetic fields is dominant, i.e. 
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where rj esc ps 0.05 between the flux of CR energy emitted 
by the shock and the flux of energy of the incoming fluid 



by the shock and t he flux ot energy or 
pu) h /2 dRiquelme and Spitkovskvl 2009) 



The higher efficiency of the B JL09 model arises from the 
fact that in this model the full pressure (~ ene rgy density) 
of the CRs is driving the instability, while in the Bell (2004) 
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(33) 



(e.g. lSkillingl 1975b . where for simplicity it is assumed that 
f(x,p) depends only on one spatial coordinate x and the 
magnitude of CR momentum, p. In Eq. ( 1331 ) the momentum 
diffu sion correspon ds to the second order Fermi acceleration 
(see lLongairll2010l) . 

The diffusion coefficients in Eq. ( l33l depend on the sta- 
tistical properties of magnetic turbulence that interacts with 
the particles. Adopting the deco mposition of compressible 
MHD turbulence in modes ( m.2i Yan and Lazarianl (12002 , 
2004) identified the fast modes as the principal modes re- 
sponsible for scattering and turbulent acceleration of CRs 
in galactic environment. Later analogous conclusions were 



reach ed for the CR acceleration in clusters of galaxies (Bru netti and Lazaria 
2007h . 

The above conclusion follows from Alfven and slow modes 
being inefficient for resonance interacti on with cosmic rays 
dChandranl200dlYan and La zarian 2002). Indeed, as we dis- 
cussed in £|2.2l the two modes are anisotropic with the anisotropy 
increasing with the decrease of the scale. The resonant inter- 
action of the CRs and Alfvenic perturbations happens when 
the CR Larmor radius is of the order of the parallel scale of 
the eddy. As the eddies at scales much less than injection 
scale are very elongated, the CR samples many uncorrelated 
ed dies, which reduces the in teraction efficiency. 

Yan and Lazarianl |2002) showed that the scattering by 



fast modes, which are isotropic (CL02), dominates (see Fig.[T9l. 
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Fig. 19 Implications of turbulence for cosmic rays. Upper left: Rate of CR scattering by Alfven waves versus CR energy. The lines at the top of the 
figure are the accepted estimates obtained for Kohnogorov turbulence. The dotted curve is from|Chandran (2000). The analytical calculations are 
given by the solid line with our numerical calculations given by crosses. Upper right: The rate of CR scattering (v) by fast modes in magnetically 
dominated plasma. The rate of scattering depends on damping of the fast waves (see , which in turn depends on the ratio of gaseous to magnetic 
pressure (/} = Pgas/Pmag)- Lower left: Individual trajectories of CRs tracked by the Mont e Carlo scattering code. B is obtained through 3-D 
simul ations of MHD turbulenc e. These calculations provide estimates of CR diffusion. From Yan and Lazarian (2002, 2004). Lower right: Testing 
of the lYan and La zarian (2008) non-linear theory with numerical simulations in Beresnvak and Lazarianl l2009f) . 



However, fast modes are subject to both collisional a nd col- ally, QLT has problems in treating scattering of particles 
lisionl ess dampin J^I which was taken into account in Yan and L agarian om entum nearly perpendicular to the magnetic field 
d2004 . (seeljones et alJl973lll978l:IVolkll973[ll975llOwensll974 : 

More recent studies of cosmic ray propagation and ac- Goldstein ! 1 19761 : iFelice and KulsrudboOlh and perpendic 
celerati on that explicitly appeal to th e effect of the fas t modes ular transpo rt (see iKota and Jokipiill2000t iMatthaeus et al 



Cassano and Brunettil(Eooi):lBrunettil(l2006h : lBrunetti «ja8 Q3)zarian 



Yan and Lazarian (2008 



Yan et al. 



2008 



include 
d2007h : 

The question of the effect of Alfven modes is not triv- 
ial, however. The calculations that were discussed above use 
a quasi-linear theory (QLT), whose domain of applicabil- 
ity is limited to very strong fields with small perturbations. 
The key assumption of QLT is that the particle's orbit is 
unperturbed, significantly limits its applicability. Addition- 



19 On the basis of weak turbulence theory, (Chandran 2005) has ar- 
gued that high-frequency fast waves, which move mostly parallel to 
magnetic field, generate Alfven waves also moving mostly parallel 
to magnetic field. We expect that the scattering by these generated 
Alfven modes to be similar to the scattering by the fast modes created 
by them. Therefore we expect that the simplified approach adopted in 
lYan and Laz arian ( 2004) and the papers that followed this one to hold. 



Various non-li near theories have been proposed to im- 
Dupredl966llVolldl973l[l975[|jones et al. 



prove 



973 



the QLT (see 



(2008 



Goldstein! 1976 ). In the recent paper o flYan and Lazarian 



henceforth YL08), a nonlinear theory (NLT) based 



on the lVolkl Jl975t> suggestion was developed. The new for 
malism was applied to two major processes of CR accelera 
tion in turbulence, namely , to gyroresonanc e and to the tran 



sient t ime damping (TTD) lFisk et alJdl974l) ; lGoldstein et al. 
dl975h . 

The efficient particle-wave interactions happen when the 
condition 



03 — k\\v/l = nQ 



0,±1,2. 



(34) 
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is satisfied. In Eq. ( f34b where CO is the wave frequency, Q = 
i2o/y is the relativistic gyration frequency, jj. = cos0, 6 is 
the pitch angle of particles. TTD corresponds to « = and 
it requires compressible perturbations. The most important 
for the gyroresonance are the interactions with n = 1 . 

Contrary to QLT which assumes that the magnitude of 
the magnetic field stay constant, NLT relaxes this assump- 
tion and allows the magnetic field to change in a smooth 
way. Due to conservation of adiabatic invariant p 2 ± jB (see 



Landau and Lifshit jl975 ) the pitch angle will gr adually vary, 
resulting in resonance broadening ( Volkl 19751) . Nonlinear 
transport (NLT) formalism is based on the replacement of 
the sharp resonance between waves and particles 5(fc||V|| — 
CO ± nQ ) from QLT to the "resonance function" R„ ( YL08): 

/<Ss2>\2 
,oo i(k«v«+nn-a)t-\klv\t 2 \—J- 

R„=3i dte V s o 

Jo 



\hAv\\ 



-exp 



(ibv/J — CO + nQ) 2 
k 2 ,Av 2 



(35) 



The width of the resonance function depends on the strength 

4vi|/v_i 



of the turbulence Aji = Av\\ /v± ~ \Jb~BjB — \/Ma). For gy- 
roresonance (n = ±1,2, ...) the result depends on whether /I 
is strongly or weakly perturbed by regular field. If fi Aji, 
the result is similar to QLT, because the exponents in Eq.(l35l> 
is close to the 5-functions. For jx < Aji the result is differ- 
ent. For instance, for the case of 90° scattering /I — > and the 
resonance happens mostly at k^ res ~ £2/Av, while in QLT 

k\\,res ~ fi/v|| 

The change of the efficiency of the gyroresonance scat- 
tering and acceleration by Alfvenic modes is not sufficient to 
affect the conclusions about the inefficiency of the process 
for low energy cosmic ray scatter ing by turbulence driven 
at large injection scales. However. lYan and Lazarianl ( 120081) 
showed that TTD gets appreciably more efficient, in contrast 
with the QLT-based result. TTD due to nonlinear scattering 
can be understood as a scattering by large-scale magnetic 
compressions arising from the slow mode perturbations. 

To test th e NLT results par t icle tra cing simulations were 
performed in Beresnvak et al.l (1201 II) . Test particle simula- 
tion is the tool extensively used to study CR scattering and 
transport (e.g. Giacalone and Jokipiil 1 999t Mace et al . 2000; 



The difference between QLT and NLT has important as- 
trophysical consequences. Indeed, in some phases, such as 
hot ISM, the fast mode is strongly damped, which makes 
vital the issue of the interaction of CRs with Alfven and 
slow modes. QLT predicts the marginal interaction of these 
modes with CR. One can view the NLT result for the TTD 
as the lower limit for the turbulence interaction with CRs. 



6.3 Complexity of turbulent acceleration 
6.3.1 A network of acceleration sites 

We can view the fragmented acceleration sites as a synthesis 
of local accelerators (turbulent quasi-perpendicular shocks, 
Strong turbulence and multi island reconnection) all of them 
share the same transport characteristics i.e. particles visit 
sporadically the localized and strong electric fields were the 
particles "trapped" for a short time (called "stickiness" in 
non-equilibrium statistical mechanics) and undergo sudden 
energy jump (called in statistical mechanics "Levy flight"). 
Vlahos et al.l ( 120041) . studied the evolution of particles in- 
side a Network of Acceleration Sites (non-linear structures 
with random size) associated with (random strength) elec- 
tric fields. They use three density probabilities mentioned 
already above to analyze the transport properties of the en- 
ergetic particles. PI represented the density probability 
for the spacial transport between the "scattering centers or 
Nodes" of the accelerators. 




Qin et al. 2002). The aforementioned studies, however, used 
synthetic data for turbulent fields which is problematic due 
to a number of reason. First of all, creating synthetic tur- 
bulence data which has scale-dependent anisotropy with re- 
spect to the local magnetic field as it is required both by 
theory and simulations has not been done so far. In addi- 
tion, synthetic data has used Gaussian statistics and delta- 
correlated fields, which is hard ly appropriate for descr iption 
of strong turbulence. Therefore lBeresnvak et al.l Q201 it) used 
the magnetic field obtained as a result of numerical MHD 
simulations. The results of this testing is shown in Figure[T9l 



Fig. 20 Sketch of the basic elements of the Network of particle accel- 
erators. A particle (spiraling line) basically follows the magnetic field 
lines (solid lines), although also undergoing drifts and travel freely a 
distance until it enters the "scattering center or acceleration Node" 
(filled circle), where it is accelerated by the local "electric field E". 
After spending a time T, inside the acceleration node it move freely 
again till it meets the new "Acceleration Node". 



The P{E) is the density probability for the energy gain of 
the particle interacting with the the acceleration sites ("scat- 
tering center or the Nodes"), and finally the density proba- 
bility P(Xj) of the time the particle spend inside the "scatter- 
ing center or the Node". 
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Particle acceleration in astrophysics is a multi-scale pro- 
cess. A number of multi-scale systems starting from large 
scale "turbulent shocks" or "strong turbulence" and go down 
to the small scale structures (multiple islands and/or current 
sheets). The fundamental question which we address next is: 
Is it possible, by following the dynamics of accelerated par- 
ticle inside the environments discussed above, to reconstruct 
their transport properties? 

Transport properties of the unified accelerator. 
The key problem on the transport properties of particles ac- 
celeration is the characteristics of the orbits in a dynam- 
ical system close to equilibrium (see Fig. l2Tb ) and non- 
equilibrium systems (see Fig. [2Tb). 



<lAzz{Az,z). The probability distribution function for a par- 
ticle at time t at the position z is 




Fig. 21 (a) Random walk in dynamical systems close to equilibrium 
(normal diffusion:trajectory in the left) (b) Random walk in dynamical 
systems far from equilibrium (anomalous diffusion: trajectory on the 
right). 



In the one hand we have the well known Brownian mo- 
tion which is very well studied and analyzed (called in the 
literature as Normal Diffusion) and in the other hand we see 
orbits following a mixture of short scale random walks in- 
terrupted with large "flights" called in statistical mechanics 
"Levy flights" and repres ent in the literature a s Anomalous 
Diffusion (see the tutorial IVlahos et al.l (12008b for more de- 
tails). We present bellow the man characteristics and the dy- 
namics of the orbits and the transport equations for these 
types of orbits. 

Normal diffusion and the Fokker Planck equation. 
Deriving the the Fokker Planck equation (or Kolmogorov 
forward equation) from the characteristics of the density prob- 
abilities discussed already is very instructive, that defines 
its limitations and shows its connection with systems close 
to equilibrium. We start from a description of diffusion in 
terms of a random walk, assuming that: (i) the mean value 
of the random walk steps can be different from zero, which 
corresponds to a systematic motion of the particles, and (ii) 
we assume that both the mean and the variance can be spa- 
tially dependent, which means that the distribution of incre- 
ments depends on the spatial location, i.e. it is of the form 



P{z -Az,t- At)q Az ^ z (Az,z- Az) dAz, (36) 



which is the Chapman-Kolmogorov equation, and where now 
q Az , z (Az,z) is the probability density for being at position z 
and making a step Az in time At. The FP equation can be 
derived as follows: we expand the integrand of Eq. (f36t in a 
Taylor-series in terms of z, so that P(z,t) = J^ABdAz, with 

A=P(z,t)-d,P(z,t)At-d z P(z,t)Az+^d z 2 P(z,t)Az 2 + ...,07) 
where we have also expanded to first order in t, 

B = q AZtZ (Az,z)- d z q AZtZ (Az,z)Az+^q Az jAz,z)Az 2 + ...(38) 

(note that the Taylor expansion is with respect to the sec- 
ond argument of q Az we expand only with respect to z, not 
though with respect to Az). Keeping al l terms up to seco nd 
order in Az, we find the FP equation dVlahos et alfeOOSh 

d,P{z,t) = -d z [V(z)P(z,t)]+d 2 [D(z)P(z,t)}, (39) 

with V (z) systematic or d rift term, and D( z) = {Az 2 ) (z)/2At 
the diffusion coefficient (IGardinerlll994l) . The FP equation 
is also applied to velocity space, e.g. in plasma physics in 
order to treat collisional effects, or to position and velocity 
space together. It has the advantage of being a deterministic 
differential equations that allows to describe the evolution 
of stochastic systems, as long as the diffusivities and drift 
velocities are known, and as long as the conditions for its 
applicability are met. 

From its derivation it is clear that the FP equation is 
suited only for systems close to equilibrium, with just small 
deviations of some particles from equilibrium, or, in the ran- 
dom walk sense, with just small steps of the particles per- 
forming the random walk. A further natural generalization 
for a diffusion equation in the approach followed here would 
be not to stop the Taylor expansion in Eqs. (T37T > and ( 1381 ) at 
second order in z, but to keep all terms, which would lead to 
the so-called Kramers-Moyal expansion. More details about 
the Fokker-Plan ck equation can be found in the literature 
Gardinerl(ll994 . 

Assuming that A (z) = and D(z) = D (constant) the FP 
equation is simple 



d t Pz,t = ^P(z,t) 



(40) 



and the solution to it in infinite space is again the Gaussian 
1 



P{z,t) 



-z l /ADt 



\/47tDt 

The square mean displacement is 
<z 2 (f)>= J z 2 P(z,t)dz = 2Dt, 
which is characteristic for the normal diffusion. 



(41) 



(42) 
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Anomalous diffusion and the Fractional Diffusion Equa- 
tion 

Normal diffusion has as basic characteristic the linear scal- 
ing of the mean square displacement of the particles with 
time, (r 2 ) ~ Dt. Many different experiments though, includ- 
ing the one shown in the previous section, reveal deviations 
from normal diffusion, in that diffusion is either faster or 
slower, and which is termed anomalous diffusion. A useful 
characterization of the diffusion process is again through the 
scaling of the mean square displacement with time, where 
though now we are looking for a more general scaling of the 
form 



? (0) 



t7 



(43) 



Diffusion is then classified through the scaling index y. The 
case 7 = 1 is normal diffusion, all other cases are termed 
anomalous. The cases y > 1 form the family of super-diffusive 
processes, including the particular case 7=2, which is called 
ballistic diffusion, and the cases y < 1 are the sub-diffusive 
processes. If the trajectories of a sufficient number of parti- 
cles inside a system are known, then plotting log < r 2 > vs 
logf is an experimental way to determine the type of diffu- 
sion occurring in a given system. 

As an illustration, let us consider a particle that is mov- 
ing with constant velocity v and undergoes no collisions and 
experiences no friction forces. It then obviously holds that 
r = vt, so that (r 2 (t)) ~ t 2 . Free particles are thus super- 
diffusive in the terminology used here, which is also the 
origin of the name ballistic for the case 7=2. Accelerated 
particles would even diffuse faster. The difference between 
normal and a anomalous diffusion is also illustrated in Fig. 
|2T1 where in the case of anomalous diffusion long "flights" 
are followed by efficient "trapping" of particles in localized 
spatial regions, in contrast to the more homogeneous picture 
of normal diffusion. 

It is to note that anomalous diffusion manifests itself not 
only in the scaling of Eq. d43l with 7 7^ 1 (which experi- 
mentally may also be difficult to be measured), but also in 
'strange' and 'anomalous' phenomena such as 'uphill' dif- 
fusion, where particles or heat diffuse in the direction of 
higher concentration, or the appearance of non-Maxwellian 
distributed particle velocities very often of power-law shape, 
and it is the main interest in this proposal. 



v R t 




V R { 



Fig. 22 Cosmic rays spiral about a reconnected magnetic field line 
and bounce back at points A and B. The reconnected regions move to- 
wards each other with the reconnection velocity Vr . The advection of 
cosmic rays entrained on magnetic field lines happens at the outflow 
velocity, which is in most cases of the order of Va- Bouncing at points 
A and B happens because either of streaming instability induced by en- 
ergetic particles or magnetic turbulence in the reconnection region. In 
reality, the outflow region gets filled in by the oppositely moving tubes 
of reconnected flux which collide only to repe at on a smaller scale the 
pattern of the larger scale reconnection. From Lazarian 1 2 0051) . 



eral, as we have demonstrated for the case of the quasi- 
perpendicular shock, the strong turbulence and the multi- 
island reconnection the transport equation for the acceler - 
ated particles becomes fractional (IBian and Brown ing 2008). 

Finally, we would like to mention that while i n GS95 

model of turbulence subdiffusion is difficult to realize Yan and Lazarianl 
(2008), the superdiffusion related to the accelerated diver- 
gence of magnetic field lines in sp ace is the expected and 
confirmed property (LV99, see also Lazarian et alj|2004 ). 



7 Acceleration of cosmic rays within LV99 model of 
reconnection 

7.1 Model of first order Fermi acceleration 



Bian and Brown ing (2008) developed a simple model for 



particle acceleration based on the ideas mentioned above 
and using the Continuous Random Wal k developed a more 
gener al the fractional diffusion equation ( Metzler and Klafter 



In what follows we discuss the first order Fermi acceleration 
which arises from volume-filling reconnection^. The LV99 
presented such a model of reconnection and observations of 
the Solar magnetic field reconnectio n support the volume- 
filled idea ( Ciaravella and Ravmondll2008l) . 

Figure [22] exemplifies the simplest realization of the ac- 
celeration within the reconnection region expected within 
20001). It is demonstrated clearly (see also the review lVlahos et all V QQ model. As a particle bounces back and forth between 
20081) that depending on the form of the probabilities P(ii),Pi{E) cM^)gmg magnetic fluxes, it gains energy through the first 

we can build a transport equation which tailored on the char- — 

. ... £ . u 1 u u-i-t- • . n n n We would like to stress that Figure QT1 exemplifies only the first 

actenstics of the above probabilities e.g. assuming that Pi .P?,"^ r , a — • . . . 

~ moment of reconnection when the fluxes are just brought together. As 

are Gaussian we recover the FP equation listed above. In the rec0 nnection develops the volume of thickness A gets filled with 
the more general case were the probabilities are more gen- the reconnected 3D flux ropes moving in the opposite directions. 
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order Fermi a cceleration described in de Gouveia Dal Pino a nd LwMdami the spectrum of accelerated energetic particles for 
d2003U2005[ henceforth GL05) (see also lLazarianll2005h . 
To derive the energy spectrum of particles one can use 



the routine way of de aling with the fi rst order Fermi accel- 



eration in shocks (see Longairl ll992l) . Consider the process 
of acceleration of Mq particles with the initial energy Eq. If 
a particle gets energy /5Eq after a collision, its energy af- 
ter m collisions is ^'"Eq. At the same time if the probability 
of a particle to remain within the accelerating region is P, 
after m collisions the number of particles gets P m Mo. Thus 
\a(M/M )/\a(E/E )=iaP/iaP and 



M 
Mo 



(44) 



For the stationary state of accelerated particles the number 
M is the number of particles having energy equal or larger 
than E, as some of these particles are not lost and are accel- 
erated further. Therefore: 



N{E)dE = const x £- 1+ ( lnP / ln ^)j£ 



(45) 



To determine P and j3 consider the following process. 
The particles from the upper reconnection region see the 
lower reconnection region moving toward them with the ve- 
locity 2Vr (see Figurel22ii. If a particle from the upper region 
enters at an angle into the lower region the expected en- 
ergy gain of the particle is 8E/E = 2VrcosO /c. For isotropic 
distribution of particles their probability function is p(9) = 
2sin0cos0ii0 and therefore the average energy gain per 
crossing of the reconnection region is 



(8E/E) = — 

c Jo 



n l 2 i Vr 
2cos 2 0sin0c/0 = 4/3 — 

c 



(46) 



An acceleration cycle is when the particles return back to the 
upper reconnection region. Being in the lower reconnection 
region the particles see the upper reconnection region mov- 
ing with the speed V^. As a result, the reconnection cycle 
provides the energy increase (8E / E) cyc [ e = 8/3(Vr/c) and 



P = E/E Q =l +8/3 (V R /c) 



(47) 



Consider the case of Vdiff <C Vr. The total number of 
particles crossing the boundaries of the upper and lower fluxes 
is 2 x 1 /4(nc), where n is the number density of particles. 
With our assumption that the particles are advected out of 
the reconnection region with the magnetized plasma outflow 
the loss of the energetic particles is 2 x Vrii. Therefore the 
fraction of energetic particles lost in a cycle is VrH / [1 /4(nc)] = 
4V R /c and 



P = 1 - 4V R /c. 

Combining Eq. (05]), (03, (08]) one gets 
N(E)dE = conshE^^dE, 



(48) 



(49) 



the case whe n the back-reaction is negligible (see GL05, cf. 
lDrurvll2012l) .R 

The first order acceleration of particles entrained on the 
contracting magnetic loop can be understood from the Liou- 
ville theorem. In the process of the magnetic tubes contrac- 
tion a regular increase of the particle's energies is expected. 
The requirement for the process to proceed efficiently is to 
keep the accelerated particles within the contracting mag- 
netic loop. This introduces limitations on the particle dif- 
fusivity perpendicular to the magnetic field direction. The 
subtlety of the point above is related to the fact that while in 
the first-order Fermi acceleration in shocks magnetic com- 
pression is important, the acceleration via the LV99 recon- 
nection process is applicable even to incompressible fluids. 
Thus, unlike shocks, it is not the entire volume that shrinks 
for the acceleration, but only the volume of the magnetic flux 
tube. Thus high perpendicular diffusion of particles may de- 
couple them from the magnetic field. Indeed, it is easy to see 
that while the particles within a magnetic flux rope depicted 
in Figure [22] bounce back and forth between the converg- 
ing mirrors and get accelerated, if these particles leave the 
flux rope fast, they may start bouncing between the mag- 
netic fields of different flux ropes which may sometimes 
decrease their energy. Thus it is important that the particle 
diffusion both in the parallel and perpendicular directions 
to the magnetic field stay different. The particle anisotropy 
which arises from particles preferentially getting accelera- 
tion in terms of the parallel momentum may also be impor- 
tant. Similarly, the first order Fermi acceleration can hap- 
pen in terms of the perpendicular momentum. This is illus- 
trated in Figure [23] There the particle with a large Larmour 
radius is bouncing back and forth between converging mir- 
rors of reconnecting magnetic field systematically getting an 
increase of the perpendicular component of its momentum. 
Both processes take place in reconnection layers. 

7.2 Simulations of the acceleration of cosmic rays by 
reconnection 

The numerical studies of the cosmic ray acceleration in re- 



conne ction regions were performed in iKowal etaL ( 201 1 



2012), where to test the mechanism the data cubes obtained 
from the models of weakly stochastic magnetic reconnec- 
tion were used. For a given snapshot we obtain a full config- 
uration of the plasma flow variables (density and velocity) 
and magnetic field. We inject test particles in such an envi- 
ronment and integrate their trajectories solving the motion 
equation for relativistic charged particles 

— (ymu) = <?(E + uxB), (50) 
dt 

21 The obtained spectral index is similar to the one of Galactic cosmic 
rays. 
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Fig. 23 Particles with a large Larmor radius gyrate about the magnetic 
field shared by two reconnecting fluxes (the latter is frequently referred 
to as "guide field". As the particle interacts with converging magne- 
tized flow corresponding to the reconnecting components of magnetic 
field, the particle gets energy gain during every gyration. 



where u is the particle velocity, y= (l — m 2 /c 2 ) is the 
Lorentz factor, m and q are particle mass and electric charge, 
respectively, and c is the speed of light. 

The study of the magnetic reconnection is done using 
the magnetohydrodynamic fluid approximation, thus we do 
not specify the electric field E explicitly. Nevertheless, the 
electric field is generated either by the flow of magnetized 
plasma or by the resistivity effect and can be obtained from 
the Ohm's equation 

E = -vxB + T]j, (51) 

where v is the plasma velocity and j = V x B is the current 
density. 

In our studies we are not interested in the acceleration by 
the electric field resulting from the resistivity effects, thus 
we neglect the last term. After incorporating the Ohm's law, 
the motion equation can be rewritten as 



— (ymu) =?[(u-v)xB]. 



(52) 



In FigureO we present the time evolution of the kinetic 
energy of the particles which have their parallel and perpen- 
dicular (red and blue points, respectively) velocity compo- 
nents accelerated for three models of reconnection. The up- 
per left panel shows the energy evolution for a 2D model 
without the guide field (as in the models studied in the pre- 
vious sections). Initially, the particles pre-accelerate by in- 
creasing their perpendicular velocity component only. Later 
we observe an exponential growth of energy mostly due to 
the acceleration of the parallel component which stops af- 
ter the energy reaches values of 10 3 — 10 4 m p (where m p is 
the proton rest mass energy). Further on, particles acceler- 
ate their perpendicular component only with smaller linear 
rate in a log-log diagram. In 2.5D case, there is also an ini- 
tial slow acceleration of the perpendicular component fol- 
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Fig. 24 Kinetic energy evolution of a group of 10 4 protons in 2D mod- 
els of reconnection with a guide field fi-=0.0 and 0. 1 (top panels, re- 
spectively). In the bottom panel a fully 3D model with initial fi-=0.0 
is presented. The colors show how the parallel (red) and perpendicu- 
lar (blue) components of the particle velocities increase with time. The 
contours correspond to values 0. 1 and 0.6 of the maximum number of 
particles for the parallel and perpendicular accelerations, respectively. 
The energy is normalized by the rest proton mass energy. The back- 
ground magnetized flow with multiple curr ent sheet layers is at time 
4.0 in Alfven time units in all models. From lKowal et al . (2011). 



lowed by the exponential acceleration of the parallel veloc- 
ity component. Due to the presence of a weak guide field, 
the parallel component accelerates further to higher energies 
at a similar rate as the perpendicular one. This implies that 
the presence of a guide field removes the restriction seen in 
the 2D model without a guide field and allows the particles 
to increase their parallel velocity components as they travel 
along the guide field, in open loops rather than in confined 
2D islands. This result is reconfirmed by the 3D model in 
the bottom panel of Figure[2U where no guide field is neces- 
sary as the MHD domain in fully three-dimensional. In this 
case, we clearly see a continuous increase of both compo- 
nents, which suggests that the particle acceleration behavior 
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changes significantly when 3D effects are considered, where 
open loops replace the closed 2D reconnecting islands. 



7.3 Implications of the acceleration via reconnection 

7.3.1 Origin of the anomalous cosmic rays 

The processes of the energetic particle acceleration in the 
process of turbulent reconnection can pre-accelerate parti- 
cles to the intermediate energies helping to solve the prob- 
lem of particle injection into shocks. It can also act as the 
principal process of acceleration. Below we present the case 
where we believe that the latter takes place. 

Since the crossing of the termination shock (TS) by Voy- 
ager 1 (VI) in late 2004 and by Voyager 2 (V2) in mid 2007 
it became clear that several paradigms needed to be revised. 
Among them was the acceleration of particles. Prior to the 
encounter of the termination shock by VI the prevailing 
view was that anomalous cosmic rays (ACRs) were accel- 
erated at the TS by diffusive shock acceleration (PSA) to 



energies 1-300 MeV/nuc (e .g. iJokipii and Giaca lone 1998 



Cummings and Stondll998l) . However, with the crossing of 



the TS by VI the energy spectrum of ACR did not unroll 
to the expected source shape: a power-law at lower energies 
with a roll off at higher energies. After 2004, both the VI 
spectrum in the heliosheath and the V2 spectrum upstream 
the TS, continued to evolve toward the expected source shape. 

To explain this paradox several mode l s were proposed. 
Among them, iMcComas and Schwadron (2006) suggested 
that at a blunt shock the acceleration site for higher energy 
ACRs would be at the flanks of the TS, where the injection 
efficiency would be higher for DS A and connection times of 
the magnetic field lines to the shoc k would be l onger , allow- 
ing acceleration to higher energies. iFisk et al.1 (120061) on the 
other hand suggested that stochastic acceleration in the tur- 
bulent heliosheath would continue to accelerate ACRs and 
that the high-energy sou rce region wou ld th us be beyond the 
TS. O ther works, such as Jokipii ( 2006 ) and Florinski and Zankl 
(2006) try to expla in the deficit o f ACRs based on a dynamic 
termination shock. Jokipii ( 2006) pointed out that a shock in 
motion on time scales of the acceleration time of the ACRs, 
days to months, would c ause the spectrum to diffe r from 
the expected DSA shape. iFlorinski and Zankl (120061) calcu- 
lated the effect of Magnetic Interacting Regions (MIRs) with 
the Termination Shock on the ACR spectral shape. They 
show that there is a prolonged period of depressed inten- 
sity in mid-energies from a single MIR. Other recent works 
have included stochastic acceleration, as well as other ef- 



fects (Moraal et al .l2006UZhangl2 006: Langn er and Po tgieter 



20061 iFerreira et al.ll2007h~ It became clear after the cross- 
ing of the TS by V2 that these models would require ad- 
justments. The observations by V2 indicate for example that 
a transient did not cause the modulation shape of the V2 



spectrum at the time of its TS crossing. When both space- 
craft were in the heliosheath in late 2007, the radial gradient 
in the 13-19 MeV/nuc ions did not appear to be caused by 
a transient. The 60-74 MeV/nuc ions have no gradient, so 
no north-south or longitudinal asymmetry is observed in the 
ACR intensities at the higher ene rgies. 

In Lazarian and Ophei (2009, henceforth LO09) we pro- 
pose an alternative model, which explains the source of ACRs 
as being in the heliosheath and we appeal to magnetic recon- 
nection as a process that can accelerate particles. LO09 ex- 
plained the origin of the magnetic field reversals that induce 
magnetic reconnection in heliosheath and heliopause. 

Indeed, it is well known that magnetic field in the he- 
liosphere change polarity and induce reconnection. For in- 
stance, a^_dieSunrotates magnetic field twists into a Parker 
spiral JParkerll958[) with magn etic fields separated by a cur- 
rent sheet (see Schattenl ll97ll) . The changes of magnetic 
field are also expected due to the Solar cycle activity. 

The question now is at what part of the heliosheath we 
expect to see reversals. The structure of the magnetic field 
in the solar wind is complex. The solar magnetic field lines 
near the termination shock are azimuthal and form a spiral 
(see Figure |25T >. We expect the reconnection and the cor- 
responding energetic particle acceleration to happen at the 
heliosheath closer to the heliopause. This explains why Voy- 
agers do not see the signatures of anomalous cosmic ray ac- 
celeration as they pass the termination s hock. Appealing to 
their model of collisionless reconnection jDrake et alj (1201 Oh 
provided a similar explanation of the origin of the anoma- 
lous cosmic rays. 

7.3.2 Acceleration of cosmic rays in heliototail 



It is known that cosmic rays arrival direction has an en- 
ergy dependent large angular scale anisotropy with an am- 
plitude of order 10~ 4 — 10~ 3 . The first comprehensive obser- 
vation of this anisotropy was provided by a network of muon 
telescopes sensi tive to sub-TeV energie s and located at dif- 
ferent latitudes dNagashima et alJ 1998). More recently, an 
anisotropy was also obse rved in the multi-TeV energy range 
by the Tibet AS y arra y (lAmenomori et al ,l2006h. Super- Kamiokande 
dGuillian et alj2007b and by MILAGRO jAbdo et alj2009b . 
and the first high statistics observation in the southern hemi- 
sphere in the 10 Te V region, is being reported by IceCube 
( Abbasi et al. 2010h . The origin of the large angular scale 
anisotropy in the cosmic rays arrival direction is still un- 
known. The structure of the local interstellar magnetic field 
is likely to have an important role. However the combined 
study of the anisotropy energy and angular dependency, its 
time modulation and angular scale structure seem to suggest 
that the observation might be a combination of multiple su- 
perimposed effects, caused by phenomenologies at different 
distances from Earth. 
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Heliopause 



Fig. 25 Upper plot. Global view of the interaction of the solar wind 
with the interstellar wind. The spiral solar magnetic field (shown in 
dark dashed lines) is shown being deflected at the heliopause. The he- 
liopause itself is being deflec ted by the in terstellar magnetic field, (fig- 
ure adapted from Suess and Nerney 2006). Lower plot. A meridional 
view of the boundary sectors of the heliospheric current sheet and how 
the opposite sectors get tighter closer to the heliopause. The thickness 
of the outflow regions in the reconnection region depends on the level 
of turbulence. From LO09. 



In this context, particular interest is derived from the ob- 
servation of a broad excess of sub-TeV cosmic rays in a por- 
tion of the sky compatib le with the direction of the helio- 



spher ic tail (or heliotail) (Nagashi ma et alJ ll998: Hal l et al 



1999). The heliotail is the region of the heliosphere down- 
stream the interstellar matter wind delimited within the he- 
liopause, i.e. the bo undary that separates the solar win d and 
interstellar plasmas ( Izmodenov and Kallenb ach 2006). The 
observed excess was attributed to some unknown anisotropic 
process connected with the heliotail (thus called tail-in ex- 
cess). The gyro-radius of sub-TeV cosmic protons is less 
than about 200 AU (in a ~ 1 jiG interstellar magnetic field), 
which is approximately the size of the heliosphere and, most 
likely, smaller than the width and length of the heliotail. 
The persistence of the cosmic ray anisotropy structure in the 
multi-TeV energy range makes it challenging to link this ob- 
servation to the heliosphere. Although the unknown size and 




Fig. 26 A meridional view of the boundary sectors of the heliospheric 
current sheet and how the opposite sectors get tighter closer to the he- 
liopause and into the heliotail. The thickness of the outflow regions in 
the reconnection region depends on the level of turbulence. The length 
of the outflow regions L dep ends on the mean geom etry of magnetic 
field and turbulence. From Lazarian and Desiati (2010). 



extension of the heliotail contributes to the uncertainty on 
the energy scale at which heliospheric influence on cosmic 
rays starts to be negligible. However, we know that the ob- 
servations of multi-TeV cosmic rays anisotropy show small 
angular scale patterns superimposed to the smooth broad 
structure of the tail-in excess, which is suggestive of a lo- 
cal origin, i.e. within the heliotail. With the same technique 
used in gamma ray detection to estimate the background 
and search for sources of gamma rays, the MILAGRO col- 
laboration discovered two localized exc ess regions in th e 
cosmic rays arrival direction distribution (lAbdo et al . 2008). 
The same excess regions were reported by the ARGO-YBJ 
air shower array d Verne tto et al.1 120 09). The strongest and 
more localized of them (with an angular size of about 10°) 
coincides with the direction of the heliotail. 

The explanation of this excess related to the accelera- 
tion of e nergetic particles in reconne ction region was sug- 
gested by Lazarian and Desiati (2010). Figure l26lrepresents 
the possible structure of th e heliotail w hich arises from the 
solar magnetic field cycles (Parker 1979). The magnetic fields 
of the opposite polarities emerge as the result of 1 1 year so- 
lar dynamo cycle. As the magnetic field is carried away by 
solar wind, the reversed magnetic field regions get accumu- 
lated in the heliotail region. This is where reconnection is 
expected to occur. 

Naturally, the actual heliotail is going to be turbulent, 
which is not represented by the idealized drawing in Fig- 
ure [26] As the Alfven speed is smaller than the Solar wind 
speed, magnetic reconnection does not change the overall 
magnetic field structure. Nevertheless, the effects of turbu- 
lence are expected to be very important from the point of 
view of magnetic reconnection and the particle acceleration 
that it entails. 

The si mulations of the heliotail are extremely challeng- 
ing (see Pogorelov et al. 2009a 2009b) and have not been 
done with the sufficient resolution and extent. While we be- 
lieve that future research will provide details necessary for 
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quantitative modeling, the schematic representation of the 
heliotail structure depicted in Figure [26] is true in terms of 
major features. 

Galactic cosmic rays entering the heliotail will get re- 
accelerated, which would affect their spectrum. The esti- 
mate of the maximal energy of protons which can be ac- 
celerated through this process can be obtained through the 
usual arguments that the Larmor radius shou ld not be large r 
than the size of the magnetized region L zone ( Longairll997 ) 



; 0.5 TeV 



B 



1 iiG 



100 AU 



(53) 



which sets the limit of energies, i.e. E max , through appealing 
to the fact that protons of larger energies cannot be confined 
by magnetic fields to experience the acceleration through 
multiple bouncing back and forth between the reconnecting 
magnetic fluxes. The magnetic field estimates vary and can 
be expected of t he order of ~ 4/iG. L TOne is in the range from 
100 to 400 AU JPogorelov et alJl2009d) . Therefore cosmic 
rays with energies up to lOTeV can be re-accelerated in the 
heliotail, explaining the observation^!. 

7.3.3 Acceleration of particles by reconnection in solar 
flares and galaxy clusters 

Acceleration is known to accompany solar flares. LV99 ex- 
plains solar flares and predicts that only insubstantial por- 
tion of energy goes into heating during the reconnection. 
The rest goes into turbulence and the acceleration of ener- 
getic particles. The resulting turbulence can accelerate the 
energetic particles dPetrosian et al. 2008 ) via the second or- 
der Fermi process (see %.2\ . but the first order Fermi accel- 
eration arising from reconnection (see £17. II ) is also present 
and may dominate. The relative role of the two processes 
requires more studies. 

LV99 reconnection was discussed as the source of the 



astrophysical sources and candidates range from the birth of 
compact objects to explosions related to gamma-r ay bursts 
(GRBs), or to events in active g alaxies (AGNs) ("e.g. lMelrosd 
2009; Kote ra and Olintol 120 1 1 ) Very high energy observa- 
tions of AGNs and GRBs with the Fermi and Swift satellites 
and ground based gamma ray observatories (HESS, VERI- 
TAS and MAGIC) on the other hand, are challenging current 
theories of particle acceleration - mostly based on stochastic 
acceleration in shocks - which have to explain how parti- 
cles are accelerated to > TeV energies in regions relatively 
small compared to the fiducial scale of their sources (Sol et 
al. 2012, in prep.). 

Traditionall y discussed predominantly in the context of 
solar flares (e.g.lDrake et al. 2006, 2009; Gordovskv v et al 



2011 



2009; 



2010; Nishizuka and Shibata 2010l ; lGordovskvv a nd Bro wning] 



Zharkova et al. 20 1 1 ) , the Earth magnetotail ( Lazarian and Opher 



:otail ( LL: 

Drake et akboiol) . and the solar wind (e.g. l L azarian and Desiai 



2010), particle acceleration in magnetic reconnection sites 
is currently also being explored in relativistic astrophysical 
environments. It has been in voked in the production of ultra 
high energy cosmic rays (e.g.lde Gouyeia Dal Pino and LazarianI 
200dl200lUKotera and Olintol 120 l ib, in let-accretion disk 



systems flde Gouveia Dal 



'ino et al 



2010al2010bl:lGiannios 



2O10 Hdel Valle et al.l201 lllPing et alj2010l). and in t he gen- 



eral framework of AGNs and GRBs ( Ostrowskil2002l;lLazarian etal 



2003 



2011 



de Gouveia Dal Pino et 



alJ2010al 



Giannios 2010; Zhang and Ya: 




Uzdenskv 201 1; Uzdenskv an d McKinnevl201 ulde Gouveia Dal Pirn 



ikv 



observ ed energetic particles in Lazarian ( 2006J) and lLazarian 



d201 If) . Magnetic fields of different polarity get into contact 
as magnetized gas is being collected into the potential well 
of the galaxy cluster. In addition, superAlfvenic turbulence 
in galaxy clusters would produce magnetic field reversals 
of magnetic fields which will also dissipate through recon- 
nection. Therefore we believe that the acceleration may be 
efficient. A quantitative study of the process is presented in 
Brunetti & Lazarian (2012, in prep.). 

7.3.4 Acceleration in AGNs 

Cosmic ray acceleration in the high energy range is still not 
fully understood. The spectrum of the highest energy cosmic 
rays (UHECRs) is consistent with an origin in extragalactic 

22 An alternative expla nation based on scattering co smic rays in the 
heliotail was suggested in lDesiati and L azarian (20111) 



201 UlMcKinnev and UzdenskvEOlX 

Magnetic reconnection events like those associated to 
solar flares can be a very powerful mechanism operating on 
accretion disks (G05). In fact, the magnetic power released 
in fast reconnection flares has been found to be more than 
sufficient to accelerate relativistic plasmons and produce the 
observed radio luminosity of the nuclear jets associated both 
to microquasars and low luminous AGNs. The observed cor- 
relation b etween the radio luminosity and the black hole 
andrBflBSiflttihe se sources, which spans 10 9 orders of magnitude 
in mass dFalcke et al.ll2004l) . is naturally explained in this 
model as simply due to the magnetic reconnection activity 
at the jet launc hing region of the accretion disk corona e of 
these sources Jde Gouveia Dal Pino et al.ll2010all2010bh . A 
similar process may also explain the observed x-ray flares in 
YSOs. 



7.4 Acceleration and Gamma Ray Bursts 

In the field of Gamma Ray Bursts (GRBs), the new ob- 
servations that followed the launching of Swift and Fermi 
satellites, while solving some old problems, have raised new 
questions. The most important unknown parameter is the ra- 
tio (a) between the Poynting flux and the matter (baryonic 
+ leptoni c) flux. In the standard fireball internal shock (IS) 
scenario ( Paczvnskil 1 9861; Ishemi and Piranl lS?90h . magnetic 
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fields are assumed not to play a dynamically important role, 
i.e. cr << 1. An alternative view is that the GRB outflow 
carries a dynamically important magnetic field component, 
i.e. cr >> 1. The GRB radiation in this case would be pow- 
ered by dissipation of the magnetic field energy in the ejecta 



Dye 

(e.g.lUsovll992tlThompsonll994l;lMeszarosandReesll997 



Pirari 



1999 



2005 



Lvutikov and Blandford 2003). Several 



recent developments suggest that we should consider seri- 
ously magnetic dominated jets as a viable option for the 
source of GRBs. First, the analogy with AGNs jets in which 
it is quite certain that the inner engine cannot accelerate 
baryonic dominated jets suggest that similarly in GRBs at 
least in its inner region the jet is Poynting flux dominated. 

Recent Fermi observation of GRB 0809 16C shows that 
the bright photosphere emission associated with a putative 
fireball is missing, which raises a challenge to the traditional 
fireball IS model and suggests that the central engine likely 
launches a Poynting-flux-dominated outflow at least for this 
burst. 

In this case, magnetic energy may be sufficient to feed 
GRBs. Magnetic reconnection was sugg ested as a compo - 
nent for GRBs more than a decade ago d Thompson! I 1 994) . 
The problem lay, however, in the intrinsic difficulty of re- 
connection as it is a very slow process in ordered fields. As 
with the case for solar flares, both a slow phase of accumu- 
lation of the oppositely directed flux and a fast bursty phase 
are required for reconne ction. Essential progress was made 
by dLazarian et al.l l2003). who proposed a new GRB model 
of self-adjusted reconnection based on the findings of fast 



The perturbations build up as the mini-shells propagate 
outward (see Fig|27]r/g/!f). At a certain point, the turbulence 
reaches the critical for a fast reconnection to take place. Re- 
connection events rapidly eject outflow, which further in- 
crease the turbulence intensity. This results in a run-away 
discharge of the magnetic field energy in a reconnection/turbulence 
avalanche. This is one ICMART event, which corresponds to 
one GRB pulse. During the magnetic field energy discharge, 
the C value drops from the original value to around unity. 

A GRB is composed of several ICMART events (i.e. 
broad pulses). The peak energy E p is expected to drop from 
high to low across each pulse. The y-ray polarization de- 
gree is also expected to drop from ~ 50 — 60% to ~ a few 
% during each pulse. The magnetic field configuration at 
the end of prompt emission is largely randomized, but still 
has an ordered component. The GRB light curves should 
have two variability components, a broad (slow) component 
related to the central engine activity, and a narrow (fast) 
component associated with the relativistic magnetic turbu- 
le nce (see Figj2Tb. More detailed discussions can be found 



Zhang and Yanl (1201 II) . 



8 Discussion 

8.1 Turbulence and cosmic ray acceleration in 
astrophysical fluid 



Turbulence is the essential part of astrophysical fluid dynam- 
ics- As astrophysical turbulence is magnetized, MHD turbu- 
reconnectionin 3D turbulent magnetic fields dLazarian and Vishrftjftk presents the most important process to be accounted 



1999). As a result of the increased turbulence the reconnec- 



tion rate increases, inducing a positive feedback which re- 
sults in the explosive reconnection. Such a process provides 
an alternative explanation to gamma ray bursts (Lazarian et 
al. 2003). More recently the mode l was elaborated and con- 
nected with observational data in lZhang and Yanl d201 ll) . 

Similar to the internal shock model, the mini-shells in- 
teract internally at the radius Ris ~ r 2 cAt. Most of these 
early collisions, however, have little energy dissipation, but 
serve to distort the ordered magnetic field lines entrained 
in the ejecta. As in other astrophysical objects, the system 
is prone to turbulence because of high Reynolds and mag- 
netic Reynolds number^ Since in GRBs, the outflow is 
highly relativistic, it is impractical for the magnetic field to 
be distorted. In any case, fast reconnection can be triggered 
in weakly stochastic magnetic field (LV99). The details re- 
quire further studies especially in the regime of relativistic 
turbulence. 

23 This is a nontrivial statement given the fact that the system is 
highly collisionless. Taking into account, nevertheless, the high mag- 
netization factor, most kinetic motions are concentrated perpendicular 
to magnetic field so that it is the much suppressed perpendicular vis- 
cosity and resistivity that should be adopted. 



for. The progress in understanding MHD turbulence calls for 
the adequate representation of the turbulence in the process 
of cosmic ray propagation and acceleration. It is regretful 
that a number of outdated ad hoc ideas about MHD turbu- 
lence are still being used while describing cosmic ray inter- 
actions with turbulence. 

The decomposition of MHD turbulence into modes (Alfven, 
fast and slow) provides a way to describe the cosmic ray - 
turbulence interactions in a theory -motivated way. For in- 
stance, for the cascade driven at large scales the Alfven modes 
get marginally important for scattering of low energy cosmic 
rays. This is contrary to the claims in most of the textbooks 
where Alfven modes are described as the principal scatter- 
ing agent. On the contrary, the scattering by fast modes is 
efficient, but the damping of fast modes must be taken into 
account. This changes the treatment of acceleration com- 
pared to the accepted one. At the same time, the ranges of 
sonic and Alfven Mach number for which the decomposi- 
tion works require further studies. 

Turbulence is also important for the acceleration of cos- 
mic rays in shocks. The generation of magnetic fields both in 
the pre-shock and post-shock plasmas are important ingredi- 
ents of forming strong enough converging mirrors between 
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Fig. 27 Left panel: A cartoon picture of the ICMART model. The typical distances and o~ values of various events are marked. Right panel: An 
example of GRB light curve that shows two variability time scales. The light curve of GRB 050607 is taken from the NASA Swift GRB archival 
data web site http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/ grb_lookup.php?grb_name=050607. The predictions of decreasing gamma- 
ray polarization degree 77 and the spectral peak energy E p within individual pulses are indicatively presented. Detailed decaying functions would 
be different depending on the details of evolution of ma gnetic field config uration, o~ value, as well as balance between heating and cooling of 
electrons. The general decreasing trend is robust, (from , Zhan g a nd Ya n. 201ll) 



which cosmic rays bounce and get accelerated. In addition, 
turbulence makes reconnection of magnetic fluxes fast and 
this induces the acceleration of cosmic rays via through yet 
another process of the First order Fermi acceleration type. 

A sober note is also due. One has to accept that mag- 
netic turbulence is far from being completely understood 
phenomenon. There are different regimes of turbulence, e.g. 
turbulence in partially ionized gas presents a pronouncedly 
different behavior compared to the standard MHD turbu- 
lence picture. Similarly, turbulence in the presence of local- 
ized sources and sinks of turbulent energy, i.e. imbalanced 
turbulence, is rather different from the balanced counterpart. 
In addition, the instabilities of cosmic ray fluid and the back- 
reaction of cosmic rays on turbulence may also be important. 
The consequences of these changes on the first order Fermi 
acceleration and the second order Fermi acceleration require 
further studies. 



8.2 Models of reconnection 

Since the introduction of the LV99 model, more traditional 
approaches to reconnection have been changed considerably. 
At the time of its introduction, the models competing with 
LV99 were some modifications of a single X-point colli- 
sionless reconnection. Those models had point-wise local- 
ized reconnection region and inevitably prescribed opening 
of the reconnection region upon the scales comparable to L 
(see Figure [TTJ. Such reconnection was difficult to realize 
in astrophysical conditions in the presence of random forc- 
ing which at high probability would collapse the opening of 
the reconnection layer. Single X-point reconnection were re- 



jected in observations of solar flares bv lCiaravella and Raymond 

boosh . 

Modern models of collisionless reconnection resemble 
the original LV99 model in a number of respects. For in- 
stance, they discuss, similarly to the LV99, the volume filled 
reconnection, although one may still wonder how this vol- 
ume filling is being a chieved in the pres ence of a single re- 
connection layer (see Drake et al. 2006). While the authors 
still talk about islands produced in the reconnection, in three 
dimensions these islands a re expected to evolve into con- 
tracting 3D loops or ropes (iDaughton e t al. 2008), which is 
similar to what is depicted in Figure [22] Thus we do not 
expect to see a cardinal difference between the first order 
Fermi processes of the ac celeration described in GL03 and 
later in lDrake et al.l(12006l) . This suggests that th e back-reaction 
of the particles calculated in lDrake et al. (2006) considering 
the firehose instability may be employed as a part of the ac- 
celeration process described in GL03. 

The departure from the idea of regular reconnection and 
introduction of magnetic stochasticity is also obvious in a 

number of the recent papers appealing to the tearing mode 

instability]^] as the process of enh ancing reconnection ( Loureiro et al. 



2009- Bhattach ariee et al.ll2009ri . The 3D loops that should 



arise as a result of this process should be able to accelerate 

24 The idea of appealing to the tearing mode as a means of en - 
hancing the reconn ecti on speed can be traced back to I Strauss] l ll988h . 
IWaelbroeckl <1989t) and lShibata and Tanumal feOOlh . LV99 showed that 
the linear growth of tearing modes is insufficient to obtain fast recon- 
nection. The new attack on the problem assumes that the non-linear 
growth of the islands due to merging provides their growth rates at the 
large scales that are larger than the direct growth of the tearing modes 
at those scales. This situation when the non-linear growth is faster than 
the linear one is rather unusual and requires further investigation. 
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energetic particles via the process described in GL03. As 
tearing modes can happen in a collisional fluid, this may po- 
tentially open another channel of reconnection in such fluid. 
The limitation of this process is that the tearing mode recon- 
nection should not be too fast as this would present prob- 
lems with explaining the accumulation of the flux prior to 
the flare. At the same time the idea of tearing reconnection 
does not have the natural valve of enhancing the reconnec- 
tion speed, which is contrary to the LV model where the 
degree of reconnection is determined by the level of turbu- 
lence. Thus the periods of slow reconnection in LV99 model 
are ensured by the low level of turbulence prior to the flare. 
We believe that tearing reconnection can act to destabilize 
the laminar Sweet-Parker reconnection layer, inducing tur- 
bulence. In this sense the new ideas on tearing instability 
may be complementary to the LV99 model. 

We note, however, that in most astrophysical situations 
one has to deal with the pre-existing turbulence, which is the 
consequence of high Reynolds number of the fluid. Such tur- 
bulence may modify or suppress instabilities, including the 
tearing mode instability. We claim that it, by itself, induces 
fast reconnection. We may claim that even if the astrophys- 
ical fluid is kept initially laminar the fluid a thick outflow 
from the reconnection region caused b y tearing is expected 
to bec ome turbulent. It was argued in lLapenta and Lazarian 
(120121) that this may be the cause of the reconnecti on explo- 
sions reported recently within M HD simulations (ILapenta 
2008UBettarini and Lapentalboid) . 



9 Summary 

The results of this review can be briefly summarized as fol- 
lows: 

• Astrophysical turbulence is ubiquitous and it plays the es- 
sential role for both acceleration of cosmic rays and mag- 
netic reconnection. 

• Advances in understanding of MHD turbulence help un- 
derstanding both cosmic ray acceleration and reconnection. 
Our understanding of MHD turbulence in its different regimes 
is still incomplete, however. 

• Turbulence makes magnetic reconnection fast. The LV99 
model that describes the process has been successfully tested 
numerically. 

• Preexisting in the media density fluctuations generate mag- 
netic field while interacting with the shock precursor. This 
increases the efficiency of cosmic ray acceleration and in 
some instances can even bootstrap the process of accelera- 
tion. 

• The turbulent second order Fermi acceleration is impor- 
tant for many environments. The standard Quasi Linear The- 
ory (QLT) is not accurate to describe the process, but it can 
be improved; the improved theory corresponds to numerical 
simulations. 



• The LV99 model of reconnection entails the first order 
Fermi acceleration of energetic particles. Numerical simu- 
lations support this conclusion. 

• A number of processes, e.g. the acceleration of anomalous 
cosmic rays within the magnetosphere, acceleration of cos- 
mic rays in magnetotail as well as in galaxy clusters can be 
attributed to the acceleration caused by magnetic reconnec- 
tion. 
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